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Abstract 

We investigate under which conditions a mixed state on a finite-dimensional multipartite quan¬ 
tum system may be the unique, globally stable fixed point of frustration-free semigroup dynamics 
subject to specified quasi-locality constraints. Our central result is a linear-algebraic necessary 
and sufficient condition for a generic (full-rank) target state to be frustration-free quasi-locally 
stabilizable, along with an explicit procedure for constructing Markovian dynamics that achieve 
stabilization. If the target state is not full-rank, we establish sufficiency under an additional condi¬ 
tion, which is naturally motivated by consistency with pure-state stabilization results yet provably 
not necessary in general. Several applications are discussed, of relevance to both dissipative quan¬ 
tum engineering and information processing, and non-equilibrium quantum statistical mechanics. 
In particular, we show that a large class of graph product states (including arbitrary thermal graph 
states) as well as Gibbs states of commuting Hamiltonians are frustration-free stabilizable relative 
to natural quasi-locality constraints. Likewise, we provide explicit examples of non-commuting 
Gibbs states and non-trivially entangled mixed states that are stabilizable despite the lack of an 
underlying commuting structure, albeit scalability to arbitrary system size remains in this case an 
open question. 


1 Introduction 

1.1 Context and motivation 

Convergence of a dynamical system to a stable equilibrium point is a hallmark of dissipative, irre¬ 
versible behavior. In particular, rigorously characterizing the nature and stability of equilibrium states 
of irreversible quantum evolutions is a longstanding problem central to both the mathematical the¬ 
ory of open quantum systems and the foundations of quantum statistical mechanics In recent 

years, renewed interest in these issues has been fueled by the growing theoretical and experimental 
significance of techniques for quantum reservoir engineering [46] and dissipative quantum control [40] 
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within Quantum Information Processing (QIP). Representative applications that benefit from engi¬ 
neered dissipation include robust quantum state preparation, with implications for steady-state entan¬ 
glement [MiiniiniiMiiii], non-equilibrium topological phases of matter HU, and ground-state cooling 
[5911641148j ; as well as open-system quantum simulation 0119], steady-state dissipation-driven quantum 
computation [Sllisij, dissipative quantum gadgets and autonomous quantum error correction min], 
along with quantum-limited sensing and amplification [mill]. 

While applications are often developed by making reference to a specific physical setting, a common 
theme is the key role played by constraints, that may restrict the allowed dynamical models and the 
extent of the available manipulations. This motivates seeking a rigorous system-theoretic framework for 
characterizing controlled open-quantum system dynamics subject to given resource constraints. In this 
work, we focus on dissipative multipartite quantum systems described by time-independent quasi-local 
(QL) semigroup dynamics, capturing the fact that, in many physically relevant scenarios, both the 
coherent (Hamiltonian) and irreversible (Lindblad) contributions to the semigroup generator must act 
non-trivially only on hnite subsets of subsystems, determined for instance by spatial lattice geometry. 
The main question we address is to determine what properties an arbitrary target state of interest must 
satisfy in order to be the unique stationary (“fixed”) point for a given QL constraint, thereby making 
the state globally QL-stabilizable in principle, in an asymptotic sense. 

In previous work [55] [BO] , we have answered this question under the assumption that the target state 
is pure, providing in particular a necessary and sufficient linear-algebraic condition for the latter to be 
stabilizable without requiring Hamiltonian dynamics. Such pure states are called purely Dissipatively 
Quasi-Locally Stabilizable (DQLS). While restricting to a pure target state is both a natural and 
adequate hrst step in the context of dissipatively preparing paradigmatic entangled states of relevance 
to QIP (such as W or GHZ states), allowing for a general mixed fixed-point is crucial for a number 
of reasons. On the one hand, since mixed quantum states represent the most general possibility, this 
is a prerequisite for mathematical completeness. On the other hand, from a practical standpoint, QL 
stabilization of a mixed state which is sufficiently close to an “unreachable” pure target may still be 
valuable for QIP purposes, a notable example being provided by thermal graph states at sufficiently 
low temperature [54] . Furthermore, as physical systems in thermal equilibrium are typically far from 
pure, characterizing mixed-state QL stabilization might offer insight into thermalization dynamics as 
occurring in Nature and on a quantum computer [48] . From this point of view, a stability analysis of 
thermal states of QL Hamiltonians is directly relevant to developing efficient simulation and sampling 
algorithms for the quantum canonical ensemble, so-called “quantum Gibbs samplers,” as analyzed in 
[29] for commuting Hamiltonians. 

In the mixed-state scenario, the problem of QL stabilization involves qualitatively different features 
and is substantially more complex. This is largely due to the fact that the analysis tools used in 
the pure-state setting do not lend themselves to a formulation where the invariance property of the 
globally defined target state translates directly at the level of QL generator components. We bypass 
this difficulty by restricting to the important class of frustration-free (FF) semigroup dynamics [151129] . 
for which global invariance of a state also implies its invariance under each QL component. Physically, 
the FF property is known to hold within standard derivations of Markovian semigroup dynamics, for 
instance based on Davies’ weak coupling limit or “heat-bath” approaches generalizing classical Glauber 
dynamics EES]. 
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1.2 Summary of results 

The paper is organized of as follows. Sec. [2] is devoted to a self-contained recollection and setup of 
the required formalism. After summarizing basic facts about open quantum system dynamics and their 
associated operator ^-algebras, we pay special attention to characterizing the structure of the fixed- 
point set that semigroup dynamics support (Sec. 12.21) . and to formalizing the relevant QL constraints 
(Sec. 12.3(1 . The problem of interest, namely deciding if a given state can be stabilized by using FF 
dynamics, is formulated in Sec. 12.41 Such states will be called FFQLS. We revisit there the main line 
of reasoning and result for DQLS pure states, as this also serves to motivate the restriction to the 
class of FF dynamics. While Sec. [2] mostly reviews known results, a slight generalization of existing 
results regarding the structural characterization of fixed points for quantum semigroups is contained in 
Theorem [2 conversely, Theorem |4] identifies conditions under which an operator *-algebra may arise as 
the fixed-point set of completely positive dynamics while preserving a desired target state. 

Sections [3] and [3] are the core sections of the paper, presenting our new results on necessary and, 
respectively, sufficient conditions for FFQLS. In order to do that, we introduce in Sec. 13.11 the central 
concept of a Schmidt span, and a key result (Lemma[3]) showing how this concept allows one to construct 
the smallest eigenspace for an operator that acts non-trivially only on a factor of a tensor-product space, 
and admits a given global eigenvector. In Sec. 13.21 by specializing this result to QL completely-positive 
trace-preserving generators, we show that requiring invariance of a target state constrains the structure 
of their fixed-point set: not only must the kernel include the smallest *-operator subalgebra that contains 
the Schmidt span and is closed with respect to a modified product operation, but it also needs to be 
invariant with respect to the action of the modular group associated to the target state. Building 
on that, we establish in Theorem [7] a necessary condition that an arbitrary target state must obey in 
order to be FFQLS, under a specified locality constraint. Two implications of this result are especially 
worth highlighting: first, as we formally show in Appendix A, in the limiting case where the target 
state is pure, the FFQLS notion naturally recovers the simpler DQLS notion. Second, whereas the 
existence of a unique stationary state is typically assumed on the basis of physical properties of the 
generator [12 HU [53], our result provides an independent complementary criterion. If the target state 
is generic (full-rank), we show in Theorem [9] that the necessary condition of Theorem[7|is, in fact, also 
sujjicient. The proof relies on a new way to show that a generator is FF and is constructive, yielding 
an explicit procedure for the synthesis of a stabilizing QL FF generator. Interestingly, this generator 
may be implemented in principle by purely dissipative means (no Hamiltonian control), establishing a 
further point of contact between the FFQLS and DQLS settings. While we conjecture that Theorem [9| 
remains valid for an arbitrary (non-full-rank) target, our current proof requires an additional “support 
condition” (see Theorem 1101) which, however, need not be obeyed in general. 

In Sec. El we provide several illustrative applications of our general mathematical framework, 
motivated by QIP and quantum statistical mechanics. We begin by showing that, unlike in the pure- 
state setting, even trivially separable mixed states may fail to be FFQLS - also implying that there are 
classical correlations that cannot be generated using constrained QL resources. We go on in Sec. l5.2l to 
characterize the most general FFQLS class of mixed states on d-dimensional subsystems (qudits) that 
may be constructed from a graph, by transforming an arbitrary product state with a QL quantum circuit 
naturally built out of it - so-called graph product states. Structurally, graph states are especially simple 
in that they may be associated to a family of QL commuting “parent Hamiltonians,” which, viewed 
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in a suitable basis, are also non-interacting (single-particle). Stepping up in complexity, a large and 
important family of states also derived from commuting QL Hamiltonians are commuting Gibbs states, 
which we examine in Sec. [Ql While these states are known to be FFQLS from the analysis in [52], 
they serve as a prime example to highlight the different focus and philosophy of our system-theoretic 
approach: on the one hand, not all Gibbs states of commuting Hamiltonians need to be FFQLS, if the 
QL is taken as an input of the problem, as natural from a control-oriented perspective; on the other 
hand, for every commuting Gibbs state there exists a physically motivated QL notion (stemming, in 
particular, from the stabilizing “Davies generator” [29] 1. relative to which our necessary conditions are 
indeed seen to hold - although, we explicitly verify this statement for a specific class in one dimension 
(Proposition |6|). 

Sec. 15.41 is devoted to exploring FFQLS in settings where no manifest commuting structure may be 
leveraged. Specifically, we construct and analyze three main examples. While in the pure-state setting, 
important examples exist of target states that are QL stabilizable despite being ground states of FF 
non-commuting parent Hamiltonians - including the paradigmatic case of the spin-1 AfHeck-Kennedy- 
Lieb-Tasaki (AKLT) Hamiltonian [51[3S] - we further exhibit a FFQLS family of pure generalized Dicke 
state on qudits, for which no FF commuting parent Hamiltonian exists. It is, however, for mixed target 
states that our general methods best reveal their potential, as no other approach is available to our 
knowledge. Beyond the commuting scenario, the problem becomes substantially more complicated, both 
in terms of verifying the required FFQLS necessary conditions, and in terms of maintaining physically 
meaningful QL constraints. We find that (full-rank) Gibbs states of a transverse-field one-dimensional 
Ising Hamiltonian can be FFQLS, provided that the locality constraint becomes weaker (corresponding 
to multi-body interactions with a growing weight) as the system size grows - thereby effectively limiting 
scalability in practice. As our final example, we construct a FFQLS family of entangled mixed states on 
four qubits. Beside serving as an explicit counter-example to the support condition mentioned above, 
this additionally demonstrates how “practical stabilization” of a pure state of interest (a GHZ state 
in this case) may be possible via stabilization of a mixed state arbitrarily close to the target. Final 
remarks conclude in Sec. jb] along with a discussion of implications and open questions. 

After Appendix A, which was already mentioned, we present in Appendix B a result on fixed states 
of convex combinations of completely positive trace-preserving maps that extends a weaker statement 
used in the main text (GorollaryO. Appendix G collects technical results related to “pseudo-pure” and 
pure Dicke states discussed in Secs. Ih.ll and Fb.d.ll respectively. 

2 Preliminaries 

2.1 Notation and background 

Gonsider a finite-dimensional Hilbert space "H, dim('H) = d, and let 18('H) be the set of linear operators 
on T-L. shall denote the adjoint oi X G iB('H), with self-adjoint operators X = representing 
physical observables. The adjoint operation corresponds to the transpose conjugate when applied to 
a matrix representation of X, with the simple transpose being denoted by X"^ and the entry-wise 
conjugation by X*. To avoid confusion, we shall use I to indicate the identity operator on 
whereas I will indicate the identity map (or super-operator) from ^{TL) to itself. We shall use X = Y 
to say that X is defined as Y. 
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The convex subset ©("H) C 55("H) of trace-one, positive-semidefinite operators, called density oper¬ 
ators^ is associated to physical states. We are concerned with state changes in the Schddinger picture 
between two arbitrary points in time, say 0 and t > 0, which are described by a completely-positive 
trace-preserving (CPTP) linear map (or quantum channel) on *B('H) [I]. A map 7t is CP if and only if 
it admits an operator-sum representation, 

pit) ='T'tiPo) = ^MkPoMl, poG^in), (1) 

k 

for some {Mk} C ^{H), and is also TP if in addition = I- The operators Mk are referred to 

as (Hellwig-)Kraus operators or operation elements [35] ■ The operator-sum representation of a CPTP 
map is not unique, and new decompositions may be obtained from unitary changes of the operators 
Mk- Dual dynamicf0 with respect to the Hilbert-Schmidt inner product on 55(7^) (Heisenberg picture) 
are associated to unital CP maps T^, that is, obeying the condition T^(I) = I. 

A continuous one-parameter semigroup of CPTP maps {7)}t>o, with To = T, characterized by the 
Markov composition property TtoTs = 7t+s, for all t, s > 0, will be referred to as a Quantum Dynamical 
Semigroup (QDS) [1]. We shall denote by C the corresponding semigroup generator, 7t = with the 
corresponding dual QDS {7)^}t>o being described by the generator It is well known that C (also 
referred to as the “Liouvillian”) can be always expressed in Lindblad canonical form [24l|39l|6], that is, 
in units where fi = 1: 

p{t)=C{pit)) = -i[H, p{t)]+'^[Lkpit)Ll-]^{LlLk, pit)}^, t>0, (2) 

k 

where H = is a self-adjoint operator associated with the effective Hamiltonian (generally resulting 
from the bare system Hamiltonian plus a “Lamb shift” term), and the Lindblad (or noise) operators 
{Lk} specify the non-Hamiltonian component of the generator, resulting in non-unitary irreversible 
dynamics. Equivalently, C defines a valid QDS generator if and only if it may be expressed in the form 
(see e.g. Theorem 7.1 in [53]) 

Cip) = £ip)-iKp-\-pK'^), K = iH+^£t{I), (3) 

where £ is a CP map and the anti-Hermitian part of k identihes the Hamiltonian operator. 

We shall denote by C{H, {Lk}) the QDS generator associated to Hamiltonian H and noise operators 
{Lk}. Throughout this paper, both H and all the Lk will be assumed to be time-independent, with ([2|) 
thus defining a linear time-invariant dynamical system. It is important to recall that, as for CPTP maps, 
the Lindblad representation is also not unique, namely, the same generator can be associated to different 
Hamiltonian and noise operators (see e.g. Proposition 7.4 in [62]), and, further to that, the separation 
between the Hamiltonian and the noise operators is not univocally defined [SS] |62|. Specifically, the 
Liouvillian is unchanged, £(77, {Lk}) = 73(77', {7^'^}), if the new operators may be obtained as (i) linear 
combinations of 77, {Lk} and the identity, L} = Lfc-|-Cfel, 77' = 77 — (7/2) X]fc(cfeTfe~CfcL[,), with Ck € C; 
or (ii) unitary linear combinations, L{ = UkiLi, 77' = 77, with U = {uki} a unitary matrix (and the 
smaller set “padded” with zeros if needed), corresponding to a change of operator-sum representation 
for £ in Eq. ([3|). 

“^While from a probabilistic and operator-algebra viewpoint it would be more natural to consider the dynamics acting 
on the states as (pre-)dual, we follow here the standard quantum physics notation as it allows for a more direct connection 
with existing work as well as a more compact notation in our context. 
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We will denote a f-closed associative subalgebra A C IB('H) generated by a set of operators 
Xi,...,Xk as = a\g{Xi,... ,Xk}. If T = T({Mfc}) and C = C{H,{Lk}) are a CP map and a 
QDS generator, then we shall let alg{T} = alg{Mfe} and alg{£} = a\g{H, Lk}, respectively. These 
algebras are invariant with respect to the change of representation in the Kraus or, respectively, Hamil¬ 
tonian and Lindblad operators since, as remarked, equivalent representations are linearly related to 
one another. Let 0 denote the orthogonal direct sum. It is well known that any f-closed associative 
subalgebra A of 58("H) admits a block-diagonal Wedderburn decomposition |19j . namely, % may be 
decomposed in an orthogonal sum of tensor-product bipartite subspaces, possibly up to a summand: 


(0^ 


lA), 


(B) 


in such a way that 


A=( 0H(7^'^O0I 


u(B) 




©Oij, 


' ’H_r, 


(4) 


(5) 


where represents the identity operator on the factor and Or the zero operator on Hr, 

respectively. Relative to the same decomposition, the commutant A' of A in 58(H), given by A' = 
{K I \Y,X] = 0, 'iX G A}, has the dual structure 


( 0 ® ^)) ® ISiHR). 


( 6 ) 


Consider now a density operator p G 2?(H) such that supp(.A) C supp(p), where for a generic 
operator space W the support is henceforth defined as supp(IK) = UoervSupp(O). It then follows that 


Ap = Ap^ = {Y\Y = p^Xp^, X€A}C B{H) 


(7) 


is a f-closed subalgebra with respect to the modified operator product Xi OpX 2 = Xip~^X 2 , where the 
inverse is in the sense of Moore-Penrose [55] if p is not invertible. In the QIP literature, an associative 
algebra like in Eq. which may be thought of as arising from a standard associative algebra A 
upon replacing each identity factor in Eq. ([5]) with a hxed matrix in each factor has been termed 
a distorted algebra [71 [Sj. In particular, we shall call Ap a p-distorted algebra, and refer to the map 
$p(X) = piXp^ as a “distortion map”. The p-distorted algebra generated by a set of operators 
Xi,... ,Xk will be correspondingly denoted by Ap = a\gp{{Xi,..., X^}). 


2.2 Fixed points of quantum dynamical semigroups 

States that are invariant (aka stationary or “fixed”) under the dissipative dynamics of interest will play 
a central role in our analysis. Let fix(T) indicate the set of fixed points of a CP map T; when 7t = 
for t > 0, then clearly fix(7t) = ker(£). In this section, we summarize relevant results on the structure 
of fixed-point sets for CPTP maps, and slightly extend them to continuous-time QDS evolutions. 

Recall that fixed points of unital CPTP maps form a f-closed algebra: this stems from the fact that 
alg{T} = alg{T^}, along with the following result (see e.g. Theorem 6.12 in [62]): 

Lemma 1 Given a CPTP map T, the commutant alg{T}^ is contained in fix(T^). In particular, if 
there exists a positive-definite state p > 0 in fix(H), then 

aig{r}' = fix(r^). (8) 
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If 'T is CPTP and unital, its dual map always admits the identity as a fixed point of full rank. It then 
follows that fix(T) = alg{T}' [37l[8]. A similar result can be established for QDS generators (Theorem 
7.2, [Sg): 

Lemma 2 Given a QDS generator £, the commutant alg{£}' is contained in the kernel of C). In 
particular, if C{p) = 0 for some p > 0, then 

alg{>C}' = ker(/:1'). (9) 

A key result to our aim is that, in general, the set of fixed points of a QDS has the structure of a 
distorted algebra. The following characterization is known for arbitrary (non-unital) CPTP maps (see 
e.g. Corollary 6.7 in [62]): 


Theorem 1 Given a GPTP map T and a full-rank fixed point p, 

fix(T) = p^ fix(T^) pA 

Moreover, with respect to the decomposition fix(T^) = 0^ we have 

e 

where p^^^ and are full-rank density operators of appropriate dimension. 

Building on the previous results, an analogous statement can be proved for QDS dynamics: 


( 10 ) 


( 11 ) 


Theorem 2 (QDS fixed-point sets, full-rank case) Given a QDS generator C and a full-rank fixed 
point p, 

ker(£) = ker(£l)p^. (12) 

Moreover, with respect to the decomposition ker(£l) = 0 we have 


= 0 Pi^^ ® 

i 



(13) 


where 


Pi 


{A) 


and 


-(s) 


are full-rank density operators of appropriate dimension. 


Proof: In order for {e^*}t>o to be a QDS, and thus a semigroup of trace-norm contractions [T], £ must 
have spectrum in the closed left-half of the complex plane and no purely imaginary eigenvalues with 
multiplicity. It is then easy to show, by resorting to its Jordan decomposition [62] . that the following 
limit exists: 

1 f* 

7^ = lim - / e^'^dr. 

t-i-oo t Jo 

Being the limit of convex combination of CPTP maps, which form a closed convex set. Too it also CPTP. 
Furthermore, Too projects onto ker(£), namely, fix(7jx)) = ker(£), and Too has only eigenvalues 0,1 with 
simple Jordan blocks. Similarly, it follows that the unital CP map T^, = 0~ooV projects onto ker(£t). 
Using these facts along with Theorem [TJ we then have: 

ker(£) = hx(7)x:,) = p5fix(7U)p5 =p^ker(£t)pT 
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The structure of the fixed point, Eq. m, follows from Theorem [T] applied to Too ■ 


□ 


The above two theorems make it clear that, given discrete- or continuous-time CPTP dynamics 
admitting a full-rank invariant state p, the fixed-point sets fix(T) and ker(/l) are a p-distorted algebra 
with structure 

= (14) 

i 

where the states are the same for every element in fix(7^) or ker(£). In addition, since p has a 
compatible block structure [Eq. (1131) ]. it is immediate to see that fix(T) and ker(£) are invariant with 
respect to the action of the linear map Ti\{X) = p^Xp~^ for any A £ C, and in particular for the 
modular group {Mi(p} |44] . The same holds for the fixed points of the dual dynamics. 

In fact, we can show that modular invariance is also a sufficient condition for a distorted algebra to 
be a fixed-point set for a CPTP map that fixes p, as relevant to the problem of designing stabilizing 
dynamics for p. In order to do this, we need a result by Takesaki |52j . which we give in its finite¬ 
dimensional formulation (adapted from [53], Theorem 9.2): 

Theorem 3 Let A be a f-cZosed subalgebra of and p a full-rank density operator. Then the 

following are equivalent: 

(i) There exists a unital CP map such that fix(£’^) = A, (f^)^ = and £{p) = p. 

(a) A is invariant with respect to Mi, that is, for every X G Al, p^Xp~^ & A. 

These conditions are equivalent to saying that the map is a conditional expectation on A that 
preserves p. We can then prove the following: 


Theorem 4 (Existence of p-preserving dynamics) Let p be a full-rank density operator. A dis¬ 
torted algebra Ap admits a CPTP map T such that fix(T) = Ap if and only if it is invariant for 

Ml. 

2 

Proof. First, notice that, if Ap is a distorted algebra, then it is invariant for M^ \i and only if 
the “undistorted” algebra A = p~^ApP~^ is invariant for Mi. This follows from the fact that Mi 
commutes with both the distortion map and its inverse. In particular, if Ap is invariant for Mi, we 
have: 

Mi(A) = Mi{p~^App~^) = p~^Mi{Ap)p~^ c p-iApp~^ = A. 

Thus, by Theorem O a unital CP projection £1 onto A exists whose adjoint preserves p. By Theorem 
(H the CPTP dual £ = T is such that fix(T) = Ap, as desired. 

To prove the other implication, it is sufficient to notice that Eq. (HU implies that Mi leaves 
A = fix(T^) invariant, and thus 

M^{Ap) = Mi{p^App^) = p^Mi{A)p^ c p^Ap^ = Ap. 


□ 


If the dynamics admit no full-rank fixed state, we may restrict to the support of a given fixed point, 
which is an invariant subspace for the Schrodinger’s-picture evolution: 


Theorem 5 (QDS fixed-point sets, general case) Given a finite-dimensional QDS generator £, 
and a maximal-rank fixed point p with H = supp(p), let L denote the reduction of C to We then 

have 


ker(£) = (ker(£'l') © O) . 


(15) 


Proof: For any p G ker(£), the subspace = supp(p) is invariant for the dynamics [57]. Assume that 
£ = C{H, {Lk}) and let n : Tt —>■ Tt denote the partial isometry onto T-L. Define the reduced (projected) 
operators p = nplfl, H = 1177111, = llLfelll. The dynamics inside H is then determined by the 

corresponding projected Liouvillian £(77, {Lk}) [SS], and p is, by construction, a full-rank state for this 
dynamics. Hence, the fixed-point set ker(£) is the distorted algebra ker(£) = p7ker(£l)p^. 

Consider now a maximal-rank fixed point, satisfying supp(p) = % = supp(ker(£)). It then follows 
from Theorem 9 in m that "H is not only invariant but also attractive for the dynamics. This means 
that 

^lim Tr^n-‘-e^‘(po)^ = 0, Vpo G DfH). 

With H being attractive, we have that ker(£) can have support only in H, and can thus be constructed 
by appending the zero operator on so that, using Theorem [2] 

ker(£) = ker(£) © O = p^ (ker(£l) © O) pA 


□ 

In the above proof, we made the construction explicit in terms of a representation £ = £(77, {Lk}) in 
order to make it clear that the result does not hold if we consider pi ker(£l) p5, since 77 need not be 
invariant for £l. Again, it follows that p admits a block decomposition as in Eq. (US, compatible with 
that of ker(£l) on its support. 

2.3 Quasi-local semigroup dynamics 

Throughout this paper, the open quantum system of interest will consist of a finite number n of distin¬ 
guishable subsystems, defined on a tensor-product Hilbert space 

n 

= dim(77a) = da, dim(77) = 77. 

a—1 

As in |58l |60] , we shall introduce quasi-locality constraints on the system’s evolution by specifying a list 
of neighborhoods, namely, groups of subsystems on which operators are allowed to “act simultaneously”. 
Mathematically, neighborhoods {Afj} may be specified as subsets of the set of indexes labeling the 
subsystems, that is, 

A/), C {!,...,n}, j = l,...,M. 

Each neighborhood induces a bipartite tensor-product structure of 77 as 

77 = 77At, ® njj., 77At, = (g) 77a, = (g) 77a. (16) 
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Likewise, with a neighborhood structure Af = {A/}} in place, any state p £ '£>{%) uniquely determines 
a list of reduced neighborhood states {pM,}'- 

PMj = Trj^^. (p), pMi e ®CHa^), j = 1, • ■ •, Af, (17) 

where Tr-^. indicates the partial trace over Quasi-local dynamical constraints may be specified 

by requiring compatibility with the bipartitions in (1161) . in the following sense: 

Definition 1 (Neighborhood operator) An operator X £ 05(7^) is a neighborhood operator relative 
to a given neighborhood structure Af if there exists j such that the action of X is non-trivial only on 
T-LUj , that is: 

X = Xjg-. (g) , 

where Ijj'. is the identity operator on TAjj.. 

A similar definition may be given for neighborhood CPTP maps and generators. The relevant quasi- 
locality notion for QDS dynamics is then the following: 

Definition 2 (QL semigroup) A QDS generator C is Quasi-Local (QL) relative to a given neighbor¬ 
hood structure Af if it may be expressed as a sum of neighborhood generators: 

C = ^Cj, Cj = ® Tjj',. (18) 

j 

Quasi-locality of a Liouvillian is well-defined, as the structural property in Eq. Cl is defined inde¬ 
pendently of a particular representation of the generator. In terms of an explicit representation, the 
above definition is equivalent to requiring that there exists some choice C = C{H, {Lk}), such that each 
Lindblad operator Lk is a neighborhood operator and the Hamiltonian may be expressed as a sum of 
neighborhood Hamiltonians, namely: 

Lk = Lk,Mj ® I H = Hj, Hj = Hjg-. 0 . 

j 

A Hamiltonian H of the above form is called a QL Hamiltonian (often “few-body,” in the physics 
literaturej^l. Mathematically, this denomination is natural given that, for the limiting case of closed- 
system dynamics, a QL Hamiltonian so defined automatically induces a QL (Lie-)group action consistent 
with Eq. (fT51) . with Lj = i adj^^ (•) = * [Hj^ •]• 

Remark 1. The above QL notion is appropriate to describe any locality constraint that may be asso¬ 
ciated with a spatial lattice geometry and finite interactions range (e.g., spins living on the vertices of a 
graph, subject to nearest-neighbor couplings). QL semigroup dynamics have also been considered under 
less restrictive assumptions on the spatial decay of interactions and yet different QL notions may be 
potentially envisioned (e.g. based on locality in “momentum space” or relative to “error weight”). The 
present choice provides the simplest physically relevant setting that allows for a direct linear-algebraic 
analysis. We stress that, due to the freedom in the representation of the QDS generator, QL semigroup 
dynamics may still be induced by Lindblad operators that are not manifestly of neighborhood form. 
In principle, it is always possible to check the QL property by verifying whether a QDS generator C 

®In particular, the notions of neighborhood Hamiltonian and QL Hamiltonian reduce to the standard uni-local and 
local ones for non-overlapping neighborhoods, see e.g. [2]. 
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has components only in the (super-)operator subspace spanned by QL generators. While it may be 
interesting to determine more operational and efficient QL criteria in specific cases, in most practical 
scenarios (e.g. open-system simulators [15]) available Lindblad operators are typically specified in a 
preferred neighborhood form from the outset. 

2.4 Quasi-local stabilizability: Prior pure-state results and frustration-free 
semigroup dynamics 

Our main focus will be on determining conditions under which a certain state of interest, p, is guaranteed 
to be invariant and the unique asymptotically stable state for some QL dynamics. Formally, an invariant 
state p G 2)(7^) for a QDS with generator £ is said to be Globally Asymptotically Stable (GAS) if 

lim e^*(po) = P, Vpo G 2)(77). (19) 

For QDS dynamics not subject to QL constraints it is known that a state is GAS if and only if it is the 
unique fixed point |62l 155] . A definition of stabilizable states relevant to our constrained setting may 
be given as follows: 

Definition 3 (QLS state) A state p G 2)(77) is Quasi-Locally Stabilizable (QLS) relative to a neigh¬ 
borhood structure Af if there exists a QL generator £ for which p is GAS. 

Existing work has so far focused on stabilizability of a pure state, with special emphasis on steady- 
state entanglement [SS] IMl [50] • While even in this case, in general, a careful balancing of Hamiltonian 
and dissipative action is essential, a simple yet non-trivial stabilization setting arises by further requiring 
that the target can be made QLS by a generator without a Hamiltonian component, namely, by using 
dissipation alone. Given the freedom in the representation of a QDS generator, in order to formalize 
this additional constraint we introduced a standard representation for a generator £{H, {Lfe}) that fixes 
a pure state p = |'I')(4'| G as in the following result (Gorollary 1 in [60]1: 

Proposition 1 If a generator £{H,{Lf;}) makes p = |'I')(\E'| GAS, then the same generator can be 
represented in a standard form £(H, {Lk}), in such a way that = h\’^), h £ M. and = 0, 

for all k. 

In the standard representation, the target G ker(Lfc) may thus be seen as a common “dark state” 
for all the noise operators, borrowing from quantum-optics terminology. With this in mind, a pure state 
p = |4')(4'| may be defined as Dissipatively Quasi-Locally Stabilizable (DQLS) if it is QLS with H = 0 
and QL noise operators {Lfc} in standard form. Notice that such a definition implies that p is invariant 
for the dynamics relative to each neighborhood, namely, p G ker(£({Lfc})), for each k. Building on this 
QL-invariance condition allows for proving the following characterization of DQLS states |5H] : 

Theorem 6 A pure state p = |4')('I'| G 2)(77) is DQLS relative to M if and only if 

supp(p) = Pi supp(pAtfc G (20) 

k 

Remark 2. The proof of the above result includes the construction of a set of stabilizing Lindblad 
operators L^ that make p DQLS, also implying that one such operator per neighborhood always suffices. 
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It is easy to show that any rescaled version of the same operators, L'j^ = also yield a stabilizing 

QL generator C = - incorporating “model ( 7 -)robustness,” in the terminology 

of [5S] . 

However, the reasoning followed for QL stabilization of a pure state does not extend naturally to 
a general, mixed target state. The main reason is that the standard form, hence the DQLS definition 
itself, do not have a consistent analogue for mixed states. A major simplification if p is pure stems from 
the fact that it is straightforward to check for invariance, directly in terms of the generator components 
(see Proposition 1 in [SS]); for general p, we seek a definition that extends the DQLS notion, and that 
similarly allows for explicitly studying what the invariance of p means at a QL level. A natural choice 
is to restrict to the class of frustration-free dynamics. That is, in addition to the QL constraint, we 
demand that each QL term in the generator leave the state of interest invariant. Formally, we define 

m- 


Definition 4 (FF generator) A QL generator C = Frustration Free (FF) relative to a 

neighborhood structure Af = {A/)} if any invariant state p € ker(£) also satisfies neighborhood-wise 
invariance, namely, p € ker(£j) for all j. 

Beside allowing for considerable simplification, FF dynamics are of practical interest because they are, 
similar to the DQLS setting, robust to certain perturbations. As in Remark 12.41 given a QL generator 
C = define a “neighborhood-perturbed” QL generator C = ^ ^ 

then £(p) = 0 implies £j(p) = 0 for each j; therefore, Xj£j{p) = 0 for each j, and thus C'{p) = 0 as 
well. Were C not FF, then the kernel of C would not be robust against such neighborhood-perturbations 
in general. With these motivations, we introduce the notion of QL stability that we analyze for the 
remainder of this paper: 

Definition 5 (FFQLS state) A state p € 2)(7t) is Frustration-Free Quasi-Locally Stabilizable (FFQLS) 
relative to a neighborhood structure Af if it is QLS with a stabilizing generator C that is FF. 

Remarkably, studying FFQLS states will allow us to recover the results for DQLS pure states as a 
special case. In fact, a pure state is FFQLS if and only if it is DQLS. Since the proof of this claim relies 
on results that we develop in the following sections, it is presented in Appendix A. To summarize, the 
above definition consists of four distinct mathematical conditions that the generator C of the dissipative 
dynamics must obey for a given target invariant state p: 

• (QDS): £ is a generator of a CPTP continuous semigroup; 

• (QL): £ is QL, that is, £ = Lj, with Lj = £a7 0 F-jj ,; 

• (GAS): p is GAS, or equivalently ker(£) = span(p); 

• (FF): £ is FF, namely ker(£j) > span(p) for all A/) . 

The problem we are interested in is to determine necessary and sufficient conditions for a given state 
to be FFQLS and, if so, to design QL FF dynamics that achieves the task. 
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3 Frustration-Free Stabilizable States: Necessary Conditions 


In this section, we derive necessary conditions for a target state to be FFQLS. Frustration-freeness 
requires such a state to be in the kernel of each neighborhood generator. We show that, if a neighborhood 
generator is to leave a global state invariant, the size and structure of its kernel are constrained; in 
general, the kernel will be larger than the span of the reduced neighborhood state (as a vector in Hilbert- 
Schmidt space). However, if the target state is to be the unique fixed point of the QDS dynamics, then 
the intersection of all the neighborhood-generator kernels must coincide with the span of the target 
state. We shall show in Section 14.21 that this condition is also sufficient for a generic (full-rank) state 
to be FFQLS. 

3.1 Linear-algebraic tools 

Recall that, given a tensor product of two inner-product spaces V = Va 0 Vb and a vector n € H, a 
Schmidt decomposition of v is any decomposition 

n = ^ \ia^ 0 bi, 
i 

where Ui G Va, bi G Vb, Xi > 0, and {oi}, {bi} are each orthonormal sets of vectors. There are two 
instances of Schmidt decomposition which are relevant in our context, both well known within QIP |42] . 
The first is the Schmidt decomposition of a bipartite pure state \tp) G 'Ha 0Hb, namely, 

It/j) = ^ Xi\a^) 0 \b^). 

i 

The second is the so-called operator-Schmidt decomposition, whereby a bipartite operator M G B{Ha0 
Hb) = B{Ha)0B{Hb) is factorized in terms of elements in the vector spaces B{Ha) and B{Hb), relative 
to the Hilbert-Schmidt inner product. Specifically, 

M = 0H,, 

i 

where Ai G B{Ha), G B(Hb), Xi > 0, and Tr(AjAj) = Tr(H|Hj) = Sij. 

Building on the concept of Schmidt decomposition, we introduce the Schmidt span: 

Definition 6 (Schmidt span) Given a tensor product of two inner product spaces V = Va 0 Vb and 
a vector v G V with Schmidt decomposition v = 0 bi, the Schmidt span of v relative to Va is 

the subspace 

^^(n) = spanjoi gVa\v = ^ai0b^, bi G VbY (21) 

i 

Without referring to a particular tensor-product decomposition, it is possible to show that the Schmidt 
span is the image of v under partial inner product: 

S^(n) = {a GVa\cl = (I^ 0 h^)v for some b G Vb}- (22) 

One example is when Va 0 Vb is a matrix space, such as 0 (where the latter factor is 

meant as a space of row vectors). In this case, the Schmidt span of a matrix W G 0 
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relative to the first factor is simply the range of the matrix W (namely, the set of all linear 

combinations of its column vectors), namely T,a{W) = range(hh). Similarly, the Schmidt span of W 
relative to the second factor is the orthogonal complement of the kernel of W or, in other words, 
the support: Tib{W) = supp(W^). Another example is when Va®Vb is a bipartite operator space, such 
as B{'Ha) ® B{'Hb)- The Schmidt span of M G B{'Ha) ® B{'Hb) relative to B{'Ha), i-e. “on A”, is the 
operator subspace Ea(M) = {Tr^ [{Ia ® B)M \, B G B{Hb)}- 

The Schmidt span is a useful tool because conditions on how a neighborhood operator is to affect 
a global state constrains how such an operator must act on the entire operator-Schmidt span of that 
state. This intuition is formalized in the following Lemma: 

Lemma 3 (Invariance of Schmidt span) Given a vector v & Va®Vb and Ma G B{Va), if (Ma 0 
IIb)^ = Au, then (Ma ® Ib)u^ = Xv' for all v' G Sa(u) ®Vb- In particular: 

span(u) < ker(M^ 0 Is) E^(u) 0 Vb < ker(M^ 0 Is), (23) 

and 

span(u) < &x{Ma 0 Is) => Ea(u) ®Vb < hx(MA 0 Is)- (24) 


Proof: Consider the Schmidt decomposition of u, u = ® where Oi G Va, hi G Vs, 7i > 0, 

and {oi}, {bi} are each orthonormal sets of vectors. Applying the eigenvalue equation for Ma to this 
yields 

liMAtti 0 = A 7*0* 0 bi. 

i i 

Multiplying both sides by Ia 0 6], where b] is the dual vector of bj , selects out the VA-factor of the fth 
term, i.e., MAai = Aos This holds for each i and any linear combination of the o^s. By definition, the 
Schmidt span of v is Sa(u) = spanjoi}. Denoting {,5*} a basis for Vb, we may write any v' G Yia{v)®Vb 
as v' = Yfij h-ijai 0 fdj. Applying (Ma 0 Is) to this we obtain 

(Ma 0 Ib)^^' = yy P-ijMAai 0 jdj = \v'. 
b 

Thus, all elements in Ea(u) 0 Vb have eigenvalue A with respect to Ma 0 Is, as claimed. Eqs. (l23ll 
and ((Ml) follow by specializing the above result to A = 0 and A = 1, respectively. □ 

3.2 Invariance conditions for quasi-local generators 

As remarked, we require the global dynamics to be FF. This simplifies considerably the analysis, as 
global invariance of the target state is possible only if the latter is invariant for each neighborhood 
generator. Therefore, we examine the properties of a neighborhood generator that ensure the target 
state p to be in its kernel. Note that if p is factorizable relative to the neighborhood structure (i.e., 
a pure or mixed product state), p is invariant as long as each factor of p if hxed. Fach such reduced 
neighborhood state can then be made not only invariant but also attractive by a neighborhood generator, 
if the reduced states are the only elements in the kernels of the corresponding C. This automatically 
makes the global factorized state also GAS. In other words, if p is factorizable, then QL stabilizability 
is guaranteed. If the target state is non-factorizable (in particular, entangled), the above scheme need 
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not work; a non-factorizable state will have some operator Schmidt spans with dimension greater than 
one. The following Corollary, which follows from Lemma [Sj illustrates the implication of quasi-locally 
fixing a state with non-trivial operator-Schmidt spans: 

Corollary 1 Let Cj = Cj\f^ he a neighborhood Liouvillian. If p € ker(£j), then it must also be 

that 

(p) 0 BiUji-.) < ker(/:j). 

Accordingly, if each neighborhood generator Cj is to fix a non-factorizable p (as is necessary for global 
invariance with FF dynamics), then each neighborhood generator must be constructed to leave invariant, 
in general, a larger space of operators - specifically, the corresponding neighborhood operator-Schmidt 
span of p. 

However, leaving only the Schmidt spans invariant is, in general, not possible if the dynamics are 
to be CPTP, since a Schmidt span need not be a distorted algebra (as required by Theorem [5]). We 
show that, in order for p to be in the kernel of a valid QL generator, it is necessary that the dynamics 
leave certain “minimal fixed-point sets” generated by the Schmidt spans invariant as well. We give the 
following: 

Definition 7 (Minimal modular-invariant distorted algebra) Let p G 'D{'H) he a density opera¬ 
tor, and W C The minimal modular-invariant distorted algebra generated by W is the smallest 

p-distorted algebra generated by W which is invariant with respect to A4i{X) = piXp~i, where the 
inverse is in the sense of Moore-Penrose if p is not full-rank. 

In the finite-dimensional case that we consider, Pp{W) can be constructed by the following iterative 
procedure: define = algp(kF), and compute 

until = F^ = Fp{W). This particular distorted algebra is the smallest structure whose invariance 
is required if the dynamics are to be CPTP: 

Lemma 4 (Minimal fixed-point sets) Let W < B{'H) be an operator subspace containing a positive- 
semidefinite operator p such that supp(p) = supp(VF). If W < fix(T) for a CPTP map T : B{TL) —>■ 
B{'H), then 

Fp{W) < fix(r). (25) 


Proof: Given the iterative construction of Fp{W), it suffices to show that if some set W C fix(T) 
includes a density operator with supp(p) = supp(VF), then both alg^(kF) C fix(T) and Xii{W) C 
hx(T), and their support is still equal to supp(p). Since T is a CP linear map, fix(T) is closed with 
respect to linear combinations and f-adjoint. We are left to show that fix(T) is closed with respect to 
the p-modified product. Consider a partial isometry V : supp(p) H, and define the reduced map 

f : H(supp(p)) ^ H(supp(p)), f{X) = V^r{VXV^)V. 

Since p is invariant, the set of operators with support contained in supp(p) = supp(LF) is an invariant 
subspace, and thus T{VV^XVV^) = VV^T{VV^XVV^)VV^. By construction, T is CP, TP, and has 
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a full-rank fixed point p = pV. It follows from Theorems [T] and 0] that fix(T’) is a p-distorted algebra; 
hence, it is closed with respect to the modified product, as well as modular-invariant. Now, if X,Y G W 
are fixed points for T, so are X = V^XV,Y = V^YV^ for T. Since their adjoint, linear combinations 
and p-distorted products are in fix(T), we have: 

nXp-^Y) = nVV'f XVV'f p-^VV'^YVV'f) = VV'<r{VXp-^YV^)VV^ 

= Vf{Xp~^Y)V^ = VXp-^YV^ 

= VXV'<Vp-^vWYV^ = Xp-^Y. 

Hence, it must be algp(IT) < fix(T), as desired, and we still have supp(alg^(lT)) = supp(p). 

On the other hand, if X € alg^(IT), then supp(AI i (X)) G supp(IH), and we have: 

T{M i{x)) = r{vv'<pixp-^vv^) = vv'<r{vp^xp-^v^)vv'< 

p2 

= Vf{piXp-^)V'<=Vp^Xp-W< 

= p^Xp~Y 

Accordingly, Mi{W) G fix(T) and supp(AIi (A)) C supp(p) as well, as desired. □ 

3.3 Prom invariance to necessary conditions for stabilizability 

In order to apply the above lemma to our case of interest, namely, finding necessary conditions for 
FFQLS, the first step is to show that the reduced neighborhood states of p may be used to generate 
the minimal p-distorted algebra containing the Schmidt span: 

Proposition 2 Given a neighborhood Afj G Af, the support of the corresponding reduced state, p/^f^ = 
Trjj-^ (p), is equal to the support of the operator-Schmidt span E(p) • 


Proof: Since pv^ G Sa(, (p)) supp(pjv'j) < supp(Ejy^. (p)). It remains to show the opposite inclusion, 
that is, by equivalently considering the complements, that kei{p_\f-) < ker(E 7 y^ (p)). Let {ip) G ker(p_v^). 
Since pv, > 0, we then have Tr {pAfj\4’){'ip\) = Tr (p(|i/')('0| 0 I)) = 0. Let {Ei} be a positive-operator 
valued measure (POVM) on Td-jj. which is informationally complete (that is, span{i?i} = BifHjj,)). 
The POVM elements sum to I, giving Tr (p(|i/')(V'| 0 Ei)) = 0. Since each term is non-negative, 
Tr (p(|^/>)('0| 0 Ei)) = 0 for all i. Letting pi = (p(10 Ei)), we can write 0 = Tr (p(|'0)('0| 0 Ei)) = 

{'4i\pi\ip). Then, pi > 0 implies Pil'ip) = 0 for all i. Since the Ei span the operator space by 

using Eq. (1^ . we have that the corresponding pi span Yjq-.^p). Hence, \'if) G ker(Ejv^ (p)). □ 

The above Proposition, together with Lemma 0] and Lemma 01 imply the following: 

Corollary 2 If a state p is in the kernel of a neighborhood generator Cj = 0T]y^., then the minimal 

fixed-point set generated by the neighborhood Schmidt span obeys 

(p)) 0 ) < ker(4). (26) 
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Proof: Assume that p G kei{£j). By LemmaO we have < ker(£j). By Proposition 

m we also know that the support of pM^ is equal to that of £j^.{p), and hence supp(pAtj 0 %,-) = 
supp(S 7 Vj (p) 0 With this and the fact that 

PMj 0 ^ (p) 0 ^ ker(/:j), 

Lemma S] implies that 

(^A,(p) 0 < ker(£j), 

or, equivalently, {£^fj{p)) 0 B{'Hjj',) < ker(£j), as desired. □ 

Summing up the results obtained on invariance so far, and recalling that uniqueness of the equilib¬ 
rium state is necessary for GAS, we have the following necessary condition: 

Theorem 7 (Necessary condition for FFQLS) A state p is FFQLS relative to the neighborhood 
structure Af only if 

span(p) = n-Fp^/s^,(p))0S(H:^,). (27) 

3 

Proof: Let p be FFQLS relative to N. Frustration-freeness of C implies that ker(£) = p|^ ker(£j). 
From QLS, we have ker(£) = span(p). Thus, FFQLS implies 

P|ker(£j) = span(p). 

3 

Corollary [2] implies that for each neighborhood, J'pj^. 0 ) < ker(£j). Hence, 

n (P)) ® ^ n ker(£j )- 

3 3 

By construction, we also have 

span(p) < Pi (Ea^ (p)) ® 

3 

Stringing together the three relationships above, we arrive at the desired result, Eq. (EH). □ 

4 Frustration-Free Stabilizable States: Sufficient Conditions 

In this section, we move from necessary conditions for FFQL stabilization to sufficient ones, by providing 
in the process a constructive procedure to design stabilizing semigroup generators. A key step will 
be to establish a property that arbitrary (convex) sums of Liouvillians enjoy, namely the fact that, 
as long as the algebras associated with individual components of the generator are contained in the 
algebra associated to the full generator, the existence of a common full-rank fixed point suffices to prove 
frustration-freeness. Drawing on this result, we will prove that the necessary condition of the previous 
section is also sufficient in the generic case where the target state is full-rank, and then separately 
address general target states. 
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4.1 A key result on frustration-free Markovian evolntions 

Consider a QDS of the form C = where individual terms need not, at this stage, correspond to 

neighborhood generators. The following general result holds: 


Theorem 8 (Common fixed points of sums of Liouvillians) Let C = be a sum of QDS 

generators, and assume that the following conditions hold: 

(i) alg{£fc} < alg{£} for each k; 

(ii) there exists a positive definite p £ ker(£) such that p £ ker(£fe) for all k. 

Then p' is invariant under C only if it is invariant under all Ck, that is: 

p £ ker(£) => p £ ker(£fe) V k. 


Proof: By linearity of C, we clearly have that ker(£) > P|j. ker(£fe). We show that under the hypotheses, 
ker(£) < nfeker(£fc), therefore effectively implying ker(£) = nfcker(£fc). By (ii), p is a full-rank state 
in ker(£) and p £ ker(£fe) for all k. Theorem [2] implies that 

ker(£) = p^'keT{C^)p^ and ker(£fc) = p'^keT{Cl.)p^, Vfc. 

Then, by Lemma [2 we also have that 


p 2 ker(£^)p 2 = p2alg{£}'p^ and p^'ker(Cl.)p^ = p^ alglCk}'p^ ■ 


In view of condition (i), the relevant commutants satisfy 


alg{£}' < alg{£fc}', Vfc. 


The above inequality may then be used to bridge the previous equalities, yielding: 

ker(£) = p^ker(£l)p^ = p^alg{£}'p^ 

lA 

ker(£fe) = piker{Cl)p^ = p^algiCkYp^, 

for all k. From this we obtain ker(£) < p|j,ker(£fc), which completes the proof. □ 


Remark 3. We note that condition (i) above, namely alg{£fc} < alg{£}, is only ever not satisfied 
due to the presence of Hamiltonian contributions in Ck- In fact, if Ck = £fc({£j,fc}) for each k va a 
given representation, then £ = J2k^k also has a purely dissipative representation £(Ufc{^i.fe})i and 
thus alg{£fc} < alg{£}. On the other hand, suppose that Ck = £fc(iLfc, {Tj.fc}), with iLfc ^ 0 in some 
representation. This implies that £ = £(iL, with H = J2k^k- In this case, since alg(iL) 

need not contain alg(iLfe), condition (i) does not hold in general. As a trivial example, consider two 
generators associated to Hi = M and H2 = —M, with M 7^ I. Clearly, {O} = alg(iL) does not contain 
alg(M). Likewise, ii H = (C^)®^ and Z is a single-qubit Pauli operator, consider QL Hamiltonians 
Hi = ZZl and H2 = IZZ. Then alg(iLi) ^ alg(iLi -I- H2)- Intuitively, this stems from the fact that 
since noise operators enter “quadratically” (bilinearly) in the QDS, they cannot cancel each other’s 
action - unlike Hamiltonians, which by linearity may “interfere” with one another. 

Interestingly, the reasoning leading to Theorem [8] also applies to CPTP maps, with the simplification 
that, since no Hamiltonian is present, condition (i) is always satisfied. Formally: 
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Corollary 3 (Common fixed points of sums of CPTP maps) Let T = '^f.Pklk be a sum of 

CPTP maps, with pk > 0 and = V there exists a positive definite p G fix(T) such that 

p G fix(7fe) for all k, then p' is invariant under 'T only if it is invariant under all Tk, that is: 

p G fix(T) p G fix(7fc). 

The proof is a straightforward adaptation of the one above, and it can actually be extended to a positive 
semi-definite fixed state p, provided that some extra hypotheses on the support of p are satisfied - see 
Appendix B for a precise statement and its proof. 

Another direct corollary of Theorem [SJ which now specializes to locality-constrained dynamics, pro¬ 
vides us with a useful tool to ensure that a QL generator be FF: the generator itself and all of its QL 
components must share a full-rank fixed state. 

Corollary 4 (Prustration-freeness from full-rank fixed point) Let C = be a QL generator, 

and assume that the following conditions hold: 

(i) 3\g{Cj} < alg{£} for each j; 

(ii) there exists a positive-definite p G ker(£) such that p G ker(£j) for all j. 

Then the QL generator C is FF. 

4.2 Sufficient conditions for fnll-rank target states 

Theorem 9 (Sufficient condition for full-rank FFQLS) A full-rank state p is FFQLS relative to 
the neighborhood structure Af if 

span(p) = (28) 

j 

Proof: To show that this condition suffices for FFQLS, we must show that there exists some QL FF 
Liouvillian C for which span(p) = ker(£). Our strategy is to first construct a QL generator for which 
p is the unique state in the intersection of the QL-components’ kernels. Then, we use Thm. 2] to show 
that this generator is FF, yielding the desired equality. 

Fix an arbitrary neighborhood Afj G Af, with associated bipartition % = 'Hjg-- ® We shall 

construct a neighborhood CPTP map £j = Sjg-^ where £j^f. projects onto the minimal fixed-point 

set containing the neighborhood-Schmidt span (that is, such projection maps are duals to a conditional 
expectation). Since, by construction, (E^y^ (p)) is a modular-invariant distorted subalgebra of 
bSifHjg-^), Theorem m ensures that there exists a CPTP map such that 

(P))- 

In particular, we take £‘jj-. = £j\fj-, so that it projects onto its hxed points. Explicitly, its structure 
follows from the decomposition in Eq. m- 

^P;V,(SA/-,(p))=0fi(H(j))®<?, 

t 

with a corresponding Hilbert space decomposition TLj,/. = 0 and a 

( 

full-rank state on - Introducing partial isometries liip : TLip —>■ the sought-after maps £_\f^ 

can be constructed as: 

ip) = 0 i^hp^Lj) ® . (29) 

e 
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It is straightforward to verify that is CPTP. Recalling Eq. (|3]), we may then define a neighbor¬ 
hood QDS generator by taking k = (I)/2 = 1/2 and letting 

= S-Hj - l-Hi, V/. (30) 

Let now £ = Cj = ® Ijj, define the QL generator of the overall dynamics. We con¬ 

structed each Cj in such a way that 

ker(/:,) = (p)) ® Vj. 

Hence, by invoking the hypothesis (Eq. (|28|) 1. it follows that p is the unique state obeying 

span(p) = P|ker(£j) < ker(£). (31) 

3 

A priori, it is still possible that span(p) < ker(£). However, since we have chosen k = the neighbor¬ 
hood generators Cj defined in Eq. (1301) do not have any Hamiltonian contribution; recalling Remark 3, 
it follows that the algebra of the global generator contains the algebra of each neighborhood generator. 


alg{/:} > algj/lj}, Vj. 

Thus, by Corollary|Tl the generator C is FF. From C being FF, it follows in turn that 

ker(£) < P|ker(£j), 

3 

which, together with Eq. dSlD, implies span(p) = ker(£), as desired. 
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4.3 Sufficient conditions for general target states 

If the target state p is not full-rank, the necessary condition of Theorem [7] may still be shown to be 
sufficient for FFQLS if an additional condition (referred to as the “support condition” henceforth) is 
also obeyed: 

Theorem 10 (Sufficient condition for general FFQLS) An arbitrary state p is FFQLS relative 
to the neighborhood structure Af if 

span(p) = Pi (Eaa, (p)) 0 ) (32) 

3 

and 

supp(p) = p supp(pAA^ 0 Ij^.). (33) 

3 


Proof: Our strategy is to use the support condition of Eq. (1^ to reduce the non-full-rank case to 
the full-rank one. As in the proof of the previous theorem, fix an arbitrary neighborhood A//, and 
consider the maps defined in Eq. (17^ . Let P_\f^ G B{'Hjy^) denote the Hermitian projector onto 
supp(Apj^^ (EtVj )(p)), and P^. = 1^/,- —Pj^j the associated orthogonal projector. In this case, we compose 
each f/vj with the corresponding map 
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where is, like both CP and TP; 

Tr (M)) = Tr(M(PA,^ + )) = Tr(M), VM £ 

With this, consider new CPTP maps given by Ej^/. o E^., whereby it follows that new neighborhood 
generators may be constructed as 

° , Ej = 0 , (35) 

with the global evolution being driven, as before, by the QL generator £ = £j. 

Define now If to be the projector onto supp(p), and consider the positive-semidefinite function 
V{t) = 1 — Tr (IIt) , T G The derivative of V along the trajectories of the generator we just 

constructed is 

l/(r)=-^Tr(n/:,(r)). 

j 

By LaSalle-Krasowskii theorem [34], the trajectories will converge to the largest invariant set contained 
in the set of r such that the above Lyapunov function D(t-) = 0. We next show that this set must have 
support only on supp(p) = supp(p 7 Vj Since V is defined on global input operators, we first 

re-express each neighborhood generator Cj in Eq. (1351) as 

Ej = Ej o Ej — I, Ej = Ejy^ 0 , Ej = Ej^. 0 , 

where we have used the property Ej °Ej = {Ejy^ o Elf.) 0 Xjj. . Additionally, let Pj = Pjyf. denote 

the projector onto supp( {Eff. (p)) 0%P. Assume now that supp(r) ^ supp(p 7 Vfc for some 

A4 G A/”, that is, Tr(TP^) > 0. By using the explicit form of the maps Elf. given in Eq. we then 

have: 

v{t) < -Tr(n/:fc(T)) 

= -Tr{n{EkoEl)T-nT) 

= -Tr {UEk{PkrPk)) - Tr (rPit) + Tr (Hr). (36) 

Since the target state p is invariant under Ek, its support is also invariant. This implies that 

Tr {UEkiPkTPk)) > Tr (n£fc(nTn)) = Tr (Br). 

Hence, the sum of the first and the third term in Eq. (1361) is less than or equal to zero. The second 
term, on the other hand, is strictly negative. This is because: (i) we assumed that Tr(rP^) > 0; 
(ii) with n < Pfe, and Ek{Pk) having the same support of Pk by construction, it also follows that 
Tr [BEk{Pk)) > 0. We thus showed that no state r with support outside of the support of p can be in 
the attractive set for the dynamics. Hence, the dynamics asymptotically converges onto the support of 
p, which is invariant for all the Ej. By restricting to this set, the maps E^ have no effect and the same 
argument of Theorem [HI shows that the only invariant set in such a subspace is span(p), as desired. □ 

Remark 4. We note that the support condition in Eq. (l33l) is indeed a natural candidate for a sufficient 
FFQLS condition, since if p is pure, it reduces to the necessary and sufficient condition of Theorem 
ID However, it is provably not necessary in general, see Sec. 15.4.21 We conjecture that the necessary 
condition of Theorem^ is in fact sufficient for general (non-full-rank) states. However, we currently 
lack a complete proof. 
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5 Illustrative Applications 


In this section, we illustrate the general framework presented thus far through a number of examples, 
with a twofold goal in mind: to both demonstrate the applicability and usefulness of the mathematical 
tools we have developed, and to gain insight into the problem of mixed-state QL stabilization, along 
with appreciating important differences from the pure-state setting. For simplicity, we shall focus 
in what follows on a multi-partite system consisting of n co-dimensional qudit subsystems, namely, 
= ®a=i^a = (C'^)®", with D = dim(H) = d". In the especially important case corresponding to 
qubit (or spin-I/2) subsystems, d = 2, we shall follow standard notation and denote by {|0), |1)} an 
orthonormal (computational) basis in €? and by {aa\a = 0,1, 2, 3} = {I, X, Y, Z} the set of single-qubit 
Pauli matrices, under the natural extension to multi-qubit operators, e.g., = Xa = I®... ... I, 

with non-trivial action occurring only on the ath factor. 

5.1 Some notable failures of quasi-local stabilizability 

Before exhibiting explicit classes of states which are provably FFQLS, it may be useful to appreciate 
some distinctive features that the mixed nature of the target state entails and, with that, the failure of 
some intuitively natural mechanisms to generate candidate FFQLS states. Recall that an arbitrary pure 
product (fully factorized) state is always DQLS (or equivalently, as shown, FFQLS) [S5] thus, in other 
words, failure of a pure target state to be FFQLS always implies some entanglement in the state. In 
contrast to that, entanglement is not necessary for failures of FFQLS if the target is mixed. Consider n 
qubits arranged on a line, with neighborhood specified by nearest-neighbor (NN) pairs, Nj = {j, j -I- I}, 
j = I,...,n—1, and the manifestly separable target state 

P^p,ep = i(|0)(0|®" + |l)(ir). (37) 

Because p is already in Schmidt decomposition form for all n, it easily follows that each Schmidt span 
has the form 

Sag (p) = span{|00)(00|yj+i, |ll)(ll|jj+i}, 

and is already a pTV'j-distorted algebra invariant for AI i. Taking the intersection over all neighborhoods 
then leaves the two-dimensional space. 

Pi Ejv, (p) 0 = span{|0)(0|®", |I)(I|®”} > span(p), 

j 

which violates the necessary condition for FFQLS of Theorem [T] 

Likewise, mixing a pure FFQLS entangled state Itjj) with a trivially FFQLS target such as the fully 
mixed state results in a “pseudo-pure” target state of the form 

P= Ppp = (1 - e)l'</')(^l + £1/2”, (38) 

which is not FFQLS in general: an explicit example may be constructed by taking |'0) to be the 
two-excitation Dicke state on n = 4 qubits, 

Di 3 = 1(0011)) = ^(10011) -h 10101) -f 10110) -h |1001) -f |1010) -h |1100)), (39) 

V6 

which was proved to be DQLS relative to the three-body neighborhoods Ni = {1,2,3}, A /2 = {2,3,4} 
in [58] (see also Sec. 15. 4. H and Appendix C for explicit calculations). 
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5.2 Quasi-local stabilization of graph product states 

Multi-qubit pure graph states are an important resource across QIP, with applications ranging from 
measurement-based quantum computation m to stabilizer quantum error-correcting codes |50) . More 
recently, thermal graph states [331154] have been shown to both provide faithful approximations of pure 
graph states for sufficiently low temperatures and to support non-trivial multipartite bound entangle¬ 
ment over a temperature range [14] . In this section, we demonstrate that a broader class of mixed graph 
states on qudits, which we refer to as graph product states, are FFQLS. 

For the special case of qubits, both pure [351ISE] and thermal [51] graph states are known to be 
stabilizable with QL FF semigroup dynamics with respect to a natural locality notion induced by the 
graph. We recover and extend these results to d > 2 and a broader class of non-thermal graph states, 
without making reference, in the mixed-state case, to properties of the Davies QDS generator which 
is typically employed under weak-coupling-limit assumptions [HIMllSl]- The key property of graph 
product states is that they can be transformed to a product form relative to a “logical subsystem 
factorization,” following a change of basis which is effected by a sequence of commuting neighborhood 
unitary transformations (a QL quantum circuit). Commutativity of the unitaries effectively reduces the 
problem of FFQLS to one of local stabilization of product states. While for graph states this observation 
allows to directly obtain QL FF stabilizing dynamics, they nevertheless serve as a first relevant example 
of FFQLS, in preparation for cases where the tools we propose become indispensable. 

We formally define general qudit graph states following m- Let G = {V,E) be a graph, where 
vertices j GV are associated to qudits and edges (j, k) G E label their allowed pairwise interactions. A 
natural neighborhood structure is derived from G, by letting the jth neighborhood Mj comprise vertex 
j along with the subset of vertices adjacent to it. Rather than associating the graph to a particular 
state, the graph is used to construct a set of commuting unitary edge-wise operators, say {C/(j fc)}. The 
product of all such unitaries, Ug = Ho. k)&E constitutes the quantum circuit which is used to map 

an input product state into, in general, an entangled stat(|§. Each is defined as the generalized 

controlled-Z transformation m associated to a symmetric qudit Hadamard matrix iJ. Hadamard 
matrices H are defined by the conditions WH = dl, H = and \hij\ = 1, where hij = 

To each Hadamard, there exists a corresponding generalized controlled-^ gate acting on two-qudits, 
defined by C^\ij) = hij\ij). With these in place, Ui^j k) = ^(j k) ® unitary transformation 

which transforms local operators to neighborhood operators on Afj is defined by 

n cjy, 

keAfj\j 

where each operator is defined on the global Hilbert space Ji, and acting non-trivially only on the 
subsystems by which it is indexec]^. Standard pure qudit graph states may be defined as [T5] I'lpo) = 
Ug 1+)®”. Similarly, we define graph product states as 

n 

PG = Ug PprodUG^ = Ug(^^ Pj'JUgK (40) 

®In the cluster model of quantum computation, qubit graph states are constructed by applying this global unitary (as 
a sequence of commuting neighbor hood-wise actions) on an initial product state 

^ Since H is not uniquely defined, the above depend on the choice of H. For readability, our notation does not 

make this explicit. For standard qubit graph states, H is the discrete Fourier transform on C^. 


23 



where each pj is an arbitrary qudit mixed state. We note that graph product states are distinct from 
(though overlapping with) so-called graph diagonal states [25], defined as those states obtained by 
applying the circuit Uq to any state Pdiag diagonal in the eigenbasis of the Hadamard matrix. 

To construct QL Lindblad operators which stabilize a graph product state pc, we may simply 
construct the local Lindblad operators which prepare each factor pj of Eq. (HOI) in the un-rotated basis, 
and then transform these Lindblad operators with Ug- Let each factor pj be diagonalized by pj = 
L? l] where 7 ] > 0 are the ordered (with i) eigenvalues of the qudit density operator pj and 

Vj the diagonalizing unitary transformation. Stabilizing Lindblad operators may then be constructed 
as follows: 


LU. = 


^3 ^3 


i){i- 




Li+i, 




where i,j = 1,... ,n and each Lindblad operator is defined on the whole "H, but by construction acts 
non-trivially only on the neighborhood J\fj . That the resulting global dynamics C = Cj are FF 
follows from the commutativity of the neighborhood-Liouvillians Cj. 

It is interesting to note that, since any pure state Pdiag in Eq. (1401) which is diagonal in the computa¬ 
tional basis is necessarily a product, arbitrary pure qudit graph states are FFQLS. In general, however, 
since Pdiag niay be separable but not necessarily of product form, mixed graph diagonal states need not 
be FFQLS (in line with similar conclusions for trivially separable states, as discussed in Sec. IQ) . 
Remark 4: Graph Hamiltonians. Pure graph states may be equivalently defined as a special class 
of stabilizer states, by assigning to each vertex in G a stabilizer generator, taken from the generalized 
Pauli group Qn for n qudits [SU]. For qubits, for example, a graph state \'4>g) niay be seen to the the 
unique ground state of a QL graph Hamiltonian FIg that is a sum of generators of Qn of the form: 

n 

Hg^Y. Hg,, = - ^ X, (g) Zfc = -Ug^ ( ^ X,) Ug. (41) 

i=i jev keXj\j 3 

By construction, Hg is a sum of commuting terms, and may be easily seen to be FF (namely, such 
that Itpc) is also the ground state of each Hgj separately). Further to that, the last equality in Eq. 
(1411) makes it clear how the graph Hamiltonian is mapped to a (strictly) local one in the “logical 
basis”, following application of the circuit Ug- A feature that becomes evident from expressing graph- 
Hamiltonians in this form, and that is not shared by more general QL commuting Hamiltonians, is the 
large degeneracy of their eigenspaces - precisely 1/d of the global space dimensioij^. Thermal graph 
states [331 jM], relative to Hamiltonians as in Eq. dni), are a special case of graph product states, 
corresponding to each qudit being in a canonical Gibbs state, namely, pj oc exp{—fdjHGj) in Eq. (1301) 
(or, pG oc UQ{®je^^^^)UG), where /dj denotes the inverse equilibrium temperature of the jth qubit. 
Thermal qubit graph states can thus provide a scalable class of mixed multiparty-entangled states. 

The construction leading to graph product states may be generalized to arbitrary situations where 
a quantum circuit arising from commuting unitary neighborhood operators may be identified, not 
necessarily stemming from a graph. That is, say that U = n(j. k)^Ey( 3 ,k)^ with each U(^j^k), as above, 
being an edge-wise operator, with \U(j,k), = 0 for all edges (j, fc), (j', k'). Then we may define 

QL-transformed product states as resulting from the action of lA on any product input state: p = 

®For d = 2, this feature is key in enabling graph-state preparation in finite time with discrete-time dynamics designed 
via splitting-subspace approaches [3]. 
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U Pprod ■ Note that Pprod is, clearly, FFQLS in the strongest sense, relative to strictly local (single¬ 
site) neighborhoods, whereas p is FFQLS relative to the structure {Afj = Ujn(i/c)^0 (b fc)}’ which is 
imposed by the circuit. More generally, a QL commuting circuit may be used to extend neighborhoods 
of some input A/in into larger neighborhoods of some output Afout, associated to weaker QL constraints. 
FFQLS states maintain their property under this type of transformation, in the following sense: 

Proposition 3 (Circuit-transformed FFQLS) Let pi^ be FFQLS relative to A/in = {A/iny}, and 
let U = Uj, with Uj = ® and \Uj, Uj'] = 0 for all j,j'. Then the output state pout = ld{p) = 

U Pin C/t is FFQLS relative to Afout = {A/ijut,*:}, where 

A/’out.fe = Afin,k U ^ [J Afj^. (42) 

AtjnAtin,i/0 

Proof: This is easily verified by constructing FFQLS dynamics for U{p). If £in = is a FFQL 

stabilizing dynamics for pin, construct a new Liouvillian by conjugation, that is, £out = U o o 
More explicitly, 

ffout{p) = Uk Lin,i ^ n ffout,k{p), 

k i j j' k k 

where the neighborhood structure of the output Liouvillian relative to the (enlarged, in general) neigh¬ 
borhoods in Eq. (1421) follows from the commutativity of the circuit unitaries Uj, as conjugation by all but 
those unitaries constrained by Eq. has no net effect. Both the spectrum and the EE property are 
preserved as U is unitary, and pout is stabilized by £out because its kernel is W(ker(£in)) = W(span(pin)). 
□ 


Physically, the obvious way to construct a QL commuting circuit is via exponentiation of commuting 
QL Hamiltonians, namely, Uj = exp(iHj), with [Hj,Hji] = 0 for all j,j'. In fact, any QL commuting 
circuit arises in this way, in the sense that a family of QL commuting Hamiltonians may always be 
associated to U, for instance by letting Hj = —iXogUj in the basis which simultaneously diagonalizes 
all circuit unitaries. 

Remark 5: Rapid mixing. As mentioned, for both pure and thermal graph states on qubits, QL 
stabilizing dynamics have been thoroughly analyzed in the literature. In particular, rigorous upper 
bounds on the mixing time have been established, showing that such states may be efficiently prepared 
- that is, the (worst-case) convergence time scales only (poly-) logarithmically with the system size 
[ 501155 ] . Remarkably, rapid mixing has been shown to both lead to stability against QL perturbations 
of the generator [TO] and to the emergence of effective area laws [Sj. These results extend naturally to 
the broader classes of graph product states and encoded product FFQLS states considered here. 

5.3 Quasi-local stabilization of commuting Gibbs states 

In this subsection, we analyze FFQLS of another class of states derived from commuting QL Hamilto¬ 
nians. Consider a Gibbs state: 

^ j 
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where each Hj is a neighborhood-operator relative to Afj G Af. If the neighborhood Hamiltonians satisfy 
[Hj,Hji] = 0 for all j,j', pp is also called a commuting Gibbs state |29j . Characterizing QL evolutions 
that have canonical Gibbs states as their unique fixed point has both implications for elucidating aspects 
of thermalization in naturally occurring dynamics and for quantum algorithms and simulation - most 
notably, in the context of quantum generalizations of Metropolis sampling . Recent work [52] has 
shown that Gibbs states of arbitrary QL commuting Hamiltonians are FFQLS, the commutativity 
property being essential to ensure quasi-locality of either the weak-coupling (Davies) generator or the 
heat-bath QDS dynamics that dissipatively prepare them. The central result therein establishes an 
equivalence between the stabilizing dynamics being gapped and the correlations in the Gibbs state 
satisfying so-called “strong clustering”, implying rapid mixing for arbitrary one-dimensional (ID) lattice 
systems, or for arbitrary-dimensional lattice systems at high enough temperature. 

It is important to appreciate that in the derivation of such results, primitivity of the QDS generator 
is assumed from the outset, and verified, along with the QL and FF properties, by making explicit 
reference to the structure of the Davies or heat-bath generator (see respectively Lemma 9 and Theorem 
10 in |29jl. Gonversely, the QL notion is not a priori imposed as a design constraint for the dynamics, 
but again emerges from the structure of the generator itself. In this sense, our framework may be 
seen to provide a complementary approach, providing in particular a necessary condition for thermal 
dynamics to be primitive relative to a specified neighborhood structure. Let us illustrate the potential 
of our approach by focusing on the simplest setting of commuting two-body NN Hamiltonians in ID. 

Proposition 4 A full-rank state p > 0 defined on a ID lattice system is FFQLS relative to neighbor¬ 
hoods Afj = {j, J + 1} if and only if 

span(p) = ^ (^) {Sjj_|_i(p)}^ [^ {Si(p)} 0 (^) J^pj, j,^^{Sfc_fc_|_i(p)}^ . 

^ j odd ^ ^ k even ' 

Proof: The proof simply follows from noting that we can group the neighborhoods {j, j -I-1} with odd 
and even j and, since neighborhoods in the same group are not overlapping, that the intersection of 
the minimal fixed-point sets corresponding to neighborhoods in the same group corresponds to their 
product. □ 

One consequence of the above simplification is the following: 

Proposition 5 Let p > 0 be FFQLS with respect to the above ID NN neighborhood structure. If, for 
any neighborhood {1,1 -|- 1}, the minimal fixed-point set is the full algebra, that is, Tp, ,^i{Si .z+i(p)} = 
B{T-Liy+i), then p must factor as p = pi,,,i 0 pi+i...n- 

Proof: Assume that p is FFQLS and J),, ,^j{Spi+i(p)} = Then, p satisfies the intersection 

condition, which is simplified to a tensor product of two intersections 

span(p) = (int left)i,,,,_/ 0 (int right);+i_..._„, (44) 

where (int left)i,...,/ stands for 

1-2 V z i-i 

J'pj,^.+dEjj+i(p)}0H('Hi)j Pi (^pi{Si(p)}0 J'p,p,+j{Efc,fe+i(p)} 

j odd ^ k even 
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and, similarly, (hit right)i+i,..._„ is given by 


B{ni+i)® j Pi f ^p,,fc+i{Sfe,fe+i(p)} 

j odd,_7=Z+2 ^ ^ k even,fc—/ + 1 

Eq. (H^ can only be satisfied if p = pi...; 0 Pi+i...„. □ 

The above proposition captures the fact that, for a state of a ID NN-coupled chain to be FFQLS, 
there are limitations to the correlations that the state can exhibit. Similar restrictions are, generically, 
sufficiently strong to prevent Gibbs states of ID commuting NN Hamiltonians to be FFQLS relative to 
NN neighborhoods. A simple example is a ID Ising Hamiltonian: 

n 

H = J2ZjZj+u n>3, (45) 

i=i 

where periodic boundary conditions are assumed. The Schmidt span for each NN pair consists of the 
space of diagonal matrices, and is closed under generation of the distorted algebra. The intersection 
of all the Schmidt spans is thus the 2^-dimensional space of diagonal matrices, implying that, for all 
temperatures and size n, the Gibbs state is not FFQLS. In particular, no FF thermal dynamics subject 
to the NN QL constraint can stabilize (or be primitive with respect to) this state. 

The thermal dynamics of the Davies or heat-bath generators that stabilize commuting Gibbs states 
are, in fact, both FF and QL, albeit relative to a different neighborhood structure than the one solely 
determined by the system’s Hamiltonian [29] . This is most transparent in the weak-coupling derivation 
of the QDS, whereby the evolution induced by this Hamiltonian, Ut = effectively “modulates” in 

time the bare system-bath neighborhood coupling operators, in turn determining the relevant Lindblad 
operators in frequency space [1]. The net effect is that Ut acts as a QL commuting circuit, resulting in 
a neighborhood structure which is expanded with respect to the one associated to H or to the coupling 
operators alone (recall Proposition [3|) . In the specific Ising example of Eq. (1151) , Davies generators are 
QL for three-body (next-to-NN, NNN for short) neighborhoods, J\fj = {j — l,j, j -I- 1}. Similarly, one 
can generalize the idea and define an enlarged “Davies QL notion”. With respect to this QL constraint, 
commuting Gibbs states may be shown to obey our necessary and sufficient condition for FFQLS, as 
expected on physical grounds: 

Proposition 6 (FFQLS commuting Gibbs states) Gibbs states of ID NN commuting Hamiltonians 
are FFQLS relative to the Davies (NNN) neighborhood structure. 



Proof: Up to normalization and letting, for convenience, j5H = Hjj+i, 0 < (3 < oo, the Gibbs 

state of Eq. (l43ll may be written as pp = = ai 2 cr 23 •.. cTn-i.n, where the 

are pairwise-commuting, invertible matrices defined that are different from the identity only in NN 
sites. Our strategy is to first compute the minimal fixed-point set iFp.^^^{Tt 23 i{p)) and its intersection 
with (Ei 23(/9)), and then, by iterating, to show that the resulting intersection is span(p). To 
compute (S 234 (p)), we first obtain the corresponding Schmidt span. Using Eq. (|2^ . Ti 234 {p) = 
span{Tr^(/ 9 M)}, for all M G BifHi) G I 234 G ..,„). Letting T 2 = 11 ^= 5 ^ can write 

B{'Hi) ® I234 ® B{'H 5 ,...,n) = span{T2(Ai GI234 G H5 0 Ce,...,„)}, 
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where A,B,C range over all matrices acting on those sites. With this parameterization, the Schmidt 
span is simplified to 

^234:ip) = Span{Tri ((Ai 0 lY)^12)CT23CT34Tr5 ((^5 0 %)ct45)} 

= CT 23 CT 34 [E 2 (cri 2 ) 0 113 ® ^4(0-45)], 

noting that 1723 and 1734 commute with all operators in this space. To calculate (S 234 (p))j we first 
obtain the reduced state 

P 234 = 0'23l734Tr234 (cri2Cr45 ■ • ■ O’n-l.n) = 0'230'34{o'2 0 I 3 ® 74 ), 
where (72 = Tr^(cri 2 ) and <74 = It follows that ^p^^^{T, 23 a{p)) has a simple structure: 

•^P234 (^234(p)) = 0'23Cr34[.7>2(S2(cri2)) 0 I3 ® ^"^4 (S4 ((T45 ))]. 

Direct calculation verifies that J-p^^^{Y, 23 a[p)) obeys the required properties of closure under the distor¬ 
tion map ‘I>p 234 and invariance under p 234 -modular action. Similarly, we have 

■^Pi23(^123(p)) = Cri2Cr23[Il2 ® J'aa ( 23 ( 0 - 34 ))]. 

Finally, we compute the intersection of these two adjacent fixed-point sets. To highlight the necessary 
structure, we write (2fc(o-i^i+i)), where k = i or i + 1. The relevant intersection is then 

{o-i 20 - 23 [Ii 2 0 ® '^(^ 4 )]} H {a- 230 - 34 [S(Hi) 0 ® ® ^t^]}■ Factoring out the common invertible 

multiple of 0 - 23 , this intersection simplifies to 

<j23[{{<Ji2) ® ® ^(^ 4 ))} n mni) 0 0 (734J-^f)]}] 

= 723[{(7i 2) n (BiHi) 0 ni)} 0 {{Bll 0 BiUA)) n (734J-45)}. 

For each intersection, notice that one argument is contained in the other, giving 

723(712 0 <JzA^t\] = 712723734[I 123 ® 

By iterating, we find that subsequent intersections simplify to 712 ... 7j_i j [Ii^,,, j_i 0 For the 

final intersection, we may take 7 „_„_|_i = 0 l„+i, where we take the (n -|- l)-th system to be trivial. 

This leads to (E„( 7 „,„+i)) = span(I„). Thus, after taking the intersection over all fixed-point sets, 
we are left with 

P|-^p,-,j+i,j+ 2 ( 24 j+i,j+ 2 (p)) = 7 i 2 ... 7 „_i,„(Ii,...,„_i 0 span(I„)) =span(p), (46) 

which, by Theorem [HI proves that these states are FFQLS. □ 

5.4 Quasi-local stabilization beyond commuting Hamiltonians 

So far, the identification of a commuting structure has played an important role in the verification of 
the FFQLS property. Thus, it is an important question to determine the extent to which “lack of com¬ 
mutativity” may hinder FFQLS. The issue is simpler and better (albeit still only partially) understood 
for pure target states, in which case families of QL stabilizable states not stemming from a commuting 
structure have been identified for arbitrary system size and complex multi-partite entanglement pat¬ 
terns. Notably, spin-1 AKLT states in ID, which are the archetypal example of a valence-bond-solid 
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state in condensed-matter physics [5], as well as a spin-3/2 (or higher) AKLT states in 2D, which provide 
a resource for universal quantum computation nulla, are unique ground states of FF anti-ferromagnet 
Hamiltonians. As such, they are FFQLS using NN, two-body dissipative dynamics [3S]- Perhaps even 
more surprisingly, the FFQLS property still holds for long-range entangled states known as Motzkin 
states [la, which are also unique ground states of FF NN spin-1 Hamiltonians and have been proved 
to (logarithmically) violate the area law. 

In what follows, we first exhibit a family of “non-commuting” FFQLS multi-qudit generalized Dicke 
states, by also including a general result linking QL stabilizability of a pure state to its ability to be 
uniquely determined by its neighborhood marginals. Focusing then on mixed target states, we construct 
and analyze two explicit (non-scalable) examples showing that commutativity of the parent Hamiltonian 
or the generating QL circuit is, as for pure states, not necessary for FFQLS in general. 


5.4.1 Pure Dicke states on qudits 

Since the main emphasis of this paper is on mixed states, most of the technical proofs of the results 
in this section are deferred to Appendix C. Similar to qubit Dicke states from quantum optics |22j . 
qudit Dicke states may be constructed by symmetrizing an n-qudit product state in which k of the 
n subsystems are “excited” to a given single-particle state, and the remaining (n — k) are in their 
“vacuum”. We may further naturally generalize by allowing multi-level excitation. Specifically, let 
{|£)},£ = 0,...,((i—1), denote an orthonormal basis in and = {tt} the symmetric (permutation) 
group on n objects. Then each generalized Dicke state is in one-to-one correspondence with a vector 
of integers, A = {kg,... ,kd-i), where each ki specifies the occupation number (multiplicity) of each 
single-qudit state. That ii§, 

|A) = |(0...0,l...l, ..., 'd-l. .d-1 )) = - - 1 ), (47) 

fco fcl fed-1 kg ki kd-1 


where ki = N and 14- permutes the subsystems according to the permutation tt. A useful fact to our 
purpose is that generalized Dicke states admit a simple Schmidt decomposition. Consider a partition 
of the system into two groups of nA and ub subsystems, respectively. It is then easy to show that the 
Schmidt decomposition of |A) is 


|A) = y= E = 

V Uj Aa-I-Ab=A * 


Ub 

As 


(48) 


In order to specify the relevant class of FFQLS generalized Dicke states, the choice of the neighbor¬ 
hood structure is crucial. The following definition captures the required feature: 


Definition 8 A neighborhood structure {A4} is connected if for any bipartition of the subsystems, 
there is some neighborhood containing subsystems from both parts. 


Our main result is then contained in the following: 

^States of this form have been recently analyzed in [45], where additionally non-uniform superpositions are also 
considered. We maintain permutation symmetry as a hallmark our generalization of standard Dicke states. 
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Proposition 7 (FFQLS Dicke states ) Given n qudits and a connected neighborhood structure Af, 
there exists a (non-factorized) FFQLS generalized Dicke state relative to Af if d(m — 1 ) > n, where m 
is the size of the largest neighborhood in Af . 

The proof (given in Appendix C) is constructive, and yields in particular the state 

Dn,7n = 1(0^. di, F . A,. ■. , d - 1. ■. d - Ij ), r = n - (d - l)(m - 1), (49) 

m—1 m—1 r 

as a non-factorized (entangled) DQLS state for given, arbitrary system size. Nonetheless, note that the 
product of the neighborhood size and qudit dimension must be scaled accordingly. For example, if Af 
is fixed to be two-body NN (hence m = 2), then the qudit dimension itself must be at least d = n. In 
this sense, the resulting family of states is non-scalable. 

Remark 6. As a particular case of Propositions [7l we recover the fact (established in [58]) that the 
state 1(0011)) = dl 4 , 3 , also previously defined in Eq. (1^ . is FFQLS with respect to the neighborhood 
structure Afi = {1, 2, 3}, A /2 = {2,3,4}. 

Generalized Dicke states are non-trivially entangled. Their multiparty correlations have the feature 
of being uniquely determined by a proper subset of all possible marginals [45], that is, of being “uniquely 
joined” m- Remarkably, an interesting connection may be made between the extent to which an 
arbitrary pure target state is uniquely determined by the set of its neighborhood-marginals and the 
FFQLS property. This is formalized in the following: 

Proposition 8 (Unique joinability) If a pure state If) is FFQLS relative to Af, then \ fi) is uniquely 
determined by its neighborhood reduced states {pa^ } = { (l'0)('0l)}- 


Proof: By contradiction, if |'0) is not uniquely joinable, then there exists a state r, not necessarily pure, 
such that T |'i/')('0| and Tr^^(r) = pAt,, for all A/)- G Af. Clearly, supp(t) < supp(pA/j) 0 BifHjj,), for 
all j, heuce by recalling Theorem [6] it also follows that 

supp(r) < P|supp(pArJ ® = span(|'0)). 

3 

On the other hand, since by assumption supp(t) 7^ span(|'0)), the above inclusion must be strict, which 
is impossible since span(|'0)) is oue-dimensional. □ 

We now establish that the class of states of Eq. (l4^ constitute genuinely non-commuting examples 
of FFQLS. To do so, we show that, with respect to a certain class of neighborhood structures, these 
states are the unique ground states of non-commuting FF QL Hamiltonians, but they cannot be the 
unique ground states of any commuting FF QL Hamiltonians. 

First, since the ground-state space of a FF Hamiltonian is simply the intersection of the ground-state 
spaces of the individual Hamiltonian terms in the sum, replaciug each such term with the projector onto 
its excited space preserves the ground space of the global FF Hamiltonian. Without loss of geuerality, we 
may then restrict attention to QL Hamiltonians consisting of sums of neighborhood-acting projectors. 
An important consequence of this simplification is that for a given FFQL Hamiltonian, the only possible 
candidates for other FFQL Hamiltonians with the same ground-state space are sums of projectors with 
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enlarged ground-state spaces, with respect to those of the given Hamiltonian. However, the following 
Lemma demonstrates that, in the two-neighborhood case, a non-commuting FFQL Hamiltonian does 
not admit a commuting enlargement with the same ground state space: 

Lemma 5 (Commuting enlargements) Assume that n°,n 2 are two non-commuting projections 
such that 

hm (n?nO)- = n^s?^o. 

n—^oo 

//ni,n2 are enlarged commuting projectors such that 

HfeH^ = n«nfe = no, fc = i,2, (50) 

then nin 2 = HgS) with Hcg > Hgg. 

Proof: Using Eq. (1501) . we have that 

n?n^ = n?nin2n^, n^n? = n^HsHin?. 

Towards proof by contradiction, assume that Hqs = nin 2 = H^s- It follows that 

n0n° = n?n^sn°, n«n? = n«n°Gsn°. 

By using the defining property of n°,n 2 , however, the right hand-side in each of the above equalities 
simplifies to H^ Hqs n° = Hgg, for fc = 1, 2. This in turn yields 

tto — h*^ — 

111 

which contradicts the non-commuting assumption. □ 

With this Lemma in place, we now verify the genuine non-commutativity of these states. 

Proposition 9 (Non-commutativity of FFQLS Dicke states) For each (non-factorized) Dicke 
state Dn^m of Eq. there exists a neighborhood structure for which Dn^m cannot be the unique 

ground state of any commuting FF QL Flamiltonian. 

Proof: Consider the (non-factorized) Dicke state Dn,m- Correspondingly, we choose any connected 
neighborhood structure M with m-body neighborhoods, whereby there is at least one neighborhood 
(A/i, say) containing a system that is not contained in any other neighborhood. Such a neighborhood 
structure always exists, and Proposition [7] ensures that Dn^m is FFQLS relative to that. From [58] . 
any FFQLS state lip) is the unique ground state of some FF QL Hamiltonian. In particular, letting 
1^) = Dn.rm One such parent Hamiltonian is 

= = (51) 

3 3 

where is the projector onto the Schmidt span of with respect to A/}. We first show that the 

F[j do not commute with one another. From the Schmidt decomposition of Dn,rm given in Eq. (1481) . 
we have HtVi = X)a where the sum extends over all choices of m symbols from the symbols 
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in Dn^m = |A). A direct calculation the shows that the , and therefore, the Hj^ of overlapping 
neighborhoods do not commute with one another. 

In order to establish the desired result, note that enlarging the neighborhoods preserves FFQLS. 
Hence, Dn,m is also FFQLS with respect to the two-neighborhood neighborhood structure with one 
neighborhood being Mi and the other neighborhood, say, Mu, being the union of the remaining neigh¬ 
borhoods of M ■ Furthermore, by building a parent Hamiltonian out of these two projectors as in Eq. 
m, the reasoning above shows that the projectors of Dn,m with respect to Mi and Mu do not commute. 
Hence, by Lemma 0 no commuting enlargement exists for {A/i, Mu}- Assume now that a commuting 
QL parent Hamiltonian exists for the original neighborhood structure. The projectors Hj, onto the 
ground state spaces of these Hamiltonians also commute. Along with FF condition, this implies that 

|■^/')('0I = nin2 ... Hi;^! = HiHcf, 

where IVu = n 2 ... Hiyv/i • Hence, the QL Hamiltonian H = (I — Hi) -|- (I — H;/) constitutes a commuting 
enlargement for {A/i, Mu}, which we showed cannot exist. This contradiction then implies that no 
commuting enlargement can exist for M . □ 

5.4.2 Non-commuting Gibbs states 

In order to demonstrate that genuinely mixed target states may also be FFQLS despite not being obvi¬ 
ously associated to a commuting structure, specific examples may be constructed in ID by considering 
a generalization of the NN Ising Hamiltonian considered in Eq. (I45p . obtained by adding a transverse 
(magnetic) field. That iJ^: 

n—1 n 

H = -Y,ZjZ,+i-gY,X„ n>4, 5 G R+. (52) 

i=i 

In particular, we consider (full-rank) Gibbs states, constructed as in Eq. (I43L as well as variants inspired 
by non-equilibrium quench protocols, wherein an initial thermal state of a Hamiltonian with given g 
is evolved under a Hamiltonian with g' ^ g. In all cases, the relevant minimal distorted algebras 
and their intersection have been numerically constructed (in Matlab) for given QL constraints, and 
Theorem [9] used to determine FFQLS. The results are found to depend sensitively on the neighborhood 
structure: Gibbs (and generalized Gibbs) states on n = 4 qubits are found to be FFQLS for three-body 
neighborhoods (as in the corresponding commuting Ising case), however extending to n = 5,6 qubits 
requires neighborhoods to be further enlarged to allow for four-body Liouvillians. 

While a direct (numerical) verification is beyond reach, this points to the possibility that the (maxi¬ 
mal) neighborhood size will have to scale extensively as n increases, thereby preventing scalable FFQLS. 
An intuitive argument in support of this is the observation that, as the size of the “neighborhood com¬ 
plements” increase, the dimensions of the extended Schmidt spans do as well; correspondingly, the 
uniqueness of their intersection, as is required for FFQLS, becomes less likely. Despite this limitation, 
these results show the general applicability of our framework. 

^^Hamiltonians such as in Eq. ll52l are exactly solvable upon mapping to a free-fermion problem, and commuting once 
expressed in terms of appropriate quasi-particles. Indeed, thermal dynamics associated to free-fermion QDS is known to 
be hypercontractive [ 54 ], in the absence of a QL “real-space” constraint as we consider. 
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5.4.3 Entangled mixed states 

As a final application, we analyze QL stabilizability of a one-parameter family of mixed entangled 
states on n = 4 qubits. Beside illustrating the full procedure needed to check if p is FFQLS and 
to construct the stabilizing maps, this example is useful for a number of reasons: first, it reinforces 
that genuinely multipartite entangled mixed states can be FFQLS; second, it explicitly shows that 
the support condition under which we proved sufficiency for general target states in Theorem [10] is 
not necessary in general; lastly, it shows how a non-FFQLS state may still admit arbitrarily close 
(in Hilbert-Schmidt space) states that are FFQLS, that is, in control-theoretic language, it may still in 
principle allow for “practical stabilization”. The family of mixed states we analyze may be parametrized 
as follows: 

p,= (l-e)|(0011))((0011)|+e|GHZ4)(GHZ4|, ee(0,l), (53) 

with neighborhoods Ni = {1,2,3}, = {2,3,4}. Here, GHZ2 is the usual GHZ state on qubits, 

that is, |GHZ()) = (|0)®" -I- ... -|- |(i — 1)®") /'/d. As established in [58], this state is not DQLS for 
any non-trivial neighborhood structura^T as one may verify by seeing that the d-dimensional space 
span{|0, 0,..., 0),..., Kd — 1, d — 1,..., d — 1)} is contained in each extended Schmidt span, and hence 
their intersection is greater than just span(|GHZ(J)). 

Let us use the notation 123|4 to denote the partition of the index set {1, 2, 3,4} in the neighborhood 
{1,2,3} and the remaining index {4}, and similarly for 1|234: the l\2‘iA-Schmidt decomposition means 
the Schmidt decomposition with respect to such bipartition. In order to construct QL FF dynamics 
which render GAS, the first step is to compute the operator Schmidt span for each neighborhood. Two 
properties aid our analysis. First, both the Dicke and GHZ components are permutation symmetric, 
so that the analysis of the 1|234 partition carries over to that of 123|4. Second, they have compatible 
Schmidt decompositions, in the sense that we can find a single operator basis in which to Schmidt- 
decompose both of them and their mixtures. Using that 

1(0011)) = x/i72|0)|(011)) + x/i72|1)|(001)), 

IGHZ^) = v/T72|0)|000) + x/T72|1)|111), 

the desired operator Schmidt decomposition is 

Pe = i|0)(0|®[(l-e)|(011))((011)| + e|000)(000|] 

+ ^|0)(1| G [(1 - e)|(011))((001)| -f e|000)(lll|] 

+ ^|1)(0| G [(1 - e)|(001))((011)| + e|lll)(000|] 

+ i|l)(l|G[(l-e)|(001))((001)|+e|lll)(lll|]. 

Let us focus on the factors relative to subsystems {234} from each term above. To compute the minimal 
fixed-point set containing this Schmidt span, we first undo the distortion of the elements of this space 
by conjugating with respect to 

P-3I = -^(|(011 ))(( 011 )| + |( 001 ))(( 001 )|) + 4 (| 000 )( 000 | + | 111 )( 111 |), 

Vl-e ve 

GHZ states are graph states, though with respect to a star graph (i.e., a central node connected to (n— 1) surrounding 
nodes). The neighborhood structure induced by this graph is trivial, in that it consists of a single neighborhood containing 
all n qubits. Hence, we treat GHZ states as separate from graph states, considering only non-trivial QL constraints. 
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where, as noted, the above inverses are taken as the Moore-Penrose inverse (with or ^ being 
replaced by 0 in the singular cases of e = 0 or 1, respectively). Conjugation of each Schmidt basis 
element with respect to this P234 removes the e-dependence, namely, 

|( 011 ))(( 011 )| + | 000 )( 000 |, |( 011 ))(( 001 )| + | 000 )( 111 |, 

|( 001 ))(( 011 )| + | 111 )( 000 |, |( 001 ))(( 001 )| + | 111 )( 111 |. 

Via a unitary change of basis, we identify computational basis elements |000), |001), etc., with vectors 

{|000), 1(011)), |111), 1(001)), |ei), |e2), leg), |e4)}, (54) 

where |ei), [62), jea), |e4) are chosen to ensure orthonormality. This transformation reveals that the 
Schmidt-span operators share a common identity factor, as in this basis they read: 

(|6)(6|(g)i) ©o, (|6)(i|©I)© o, (|i)(o|©I) ©o, (|i)(i|©i) ©o, 

where the sector on which the zeros act is span{|ei),..., |e4)}. The span of these operators is closed 

under ^-algebra operations and constitutes a representation of the Pauli algebra. Thus, the distorted 

Schmidt spans of are already *-closed algebras. A simple calculation verifies that each distorted- 

algebra basis element is also unchanged hy Mi. To find the minimal fixed-point sets, we need to 

1 

apply the distortion map again, by conjugating the generators with PI34. We can write the Schmidt 
decomposition, with respect to the basis in Eq. (l54l) . as 

|o)(o|©r©o, |6)(i|©T©o, |i)(o|©T©o, |i)(i|©T©o, 

where we have defined r = e|0)(0| + (1 — e)|l)(l|. 

The last step is to construct QL Liouvillians for each neighborhood. The requirement is that the 
kernels of each of these are the corresponding minimal fixed-point sets. This can be obtained by 
considering the operators: 

To = |6)(i| ©I©I, T+= |6)(6| ©I©r 5 | 0 )(i|, L_ = | 6 )( 6 |©I©t 5 |i)( 0 |. 

The first Lindblad operator Lq is responsible for asymptotically preparing the subspace 
span{|000), |(011)), |111), |(001))}, while T-|_ and L_ stabilize the r factor. All three Lindblad operators 
must commute with the distorted algebra in order that it be preserved. Using the standard definitions 
of ladder operators, M = |0)(1| = (cr“)^, we rewrite the r-preparing Lindblad operators back in the 
original basis, in terms of standard Pauli matrices, as 

T+ = [cr^cr:^cr|(cr" ©CT;^ +cr4') + (T^(T;^cr4'(cr^+(Tj'+cr|)], 

L- = i/l-e [(ct^ + erg" + (74") cr^cr;^cr4^ + (ct^ + CT;^ + (74"], 

where now L_ = ^L^. Defining |(001)^) = (|001) + w|010)+a;^|100))/-\/3, uj = and similar 

terms to denote symmetric basis elements for the four dimensional space orthogonal to the symmetric 
subspace, the third Lindblad operator reads 

Lo = | 000 )(( 001 )^| + |( 011 ))(( 001 )^ 2 | + | 111 )(( 011 )<,| + |( 001 ))(( 011 )^ 2 |. 
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These Lindblad operators form the neighborhood Liouvillian on systems 234, namely, 

^234{p) = LqpLq + L^pL]^ + L_pLl_ — -{LqLq + L\L+ + p}, 

The global generator £ is obtained by constructing £123 in an analogous way, and by letting £ = 

•^234 + •^123- 

Using Matlab, we have verified that these dynamics are FF and stabilize the kernel of £ is equal 
to span(pE), as desired. In the limiting cases of e = 0,1, the above dynamics fail to have a unique fixed 
state. As we have already established, for the Dicke-state case of e = 0, FFQL stabilizing dynamics 
can be constructed by a separate procedure, whereas for the GHZ case of e = 1, no FFQL stabilizing 
dynamics exists relative to the given neighborhoods. 

Remark 7: Failure of support condition. In Theorem ITOl in order to obtain a general sufficient 
condition for FFQLS, we supplemented the necessary condition with the “support condition” of Eq. 
(PI) . We now show that pe fails the support condition despite being FFQLS. The support of pe is 
spanned by just |(00II)) and IGHZ 2 ). On the other hand, we also have 

supp(pi 23 OI 4 ) = span{|0000), 10001), |(001))|0),|(00I))|1), 

|(011))|0),|(011))|1), 11110), 11111)}, 

supp(Ii 0 P 234 ) = span{| 0000 ), 11000 ), | 0 )|( 001 )),| 1 )|( 001 )), 

10)1(011)), |1)|(011)), 10111), 11111)}. 

The intersection is found to be 

span{|0000), |(0001)), |(0011)), |(0111)), |1111)}. 

Since this intersection properly contains the support of pe, the claim follows. 

Remark 8. Since stabilization of pe is possible for e arbitrarily close to the GHZ-value of one, the 
above example demonstrates that there exist sequences of FFQLS mixed states that converge to a 
non-FFQLS pure state. While this allows for practical stabilization [34] in principle, we expect that 
the Liouvillian spectral gap will close as the non-FFQLS state is approached - making stabilization 
inefficient. If we normalize the Liouviallians to ||£||2 = 1, the corresponding gaps A are found to 
behave as A Ri 0.049 (1 — e). It remains an interesting question for future investigation to determine 
whether similar conclusions about QL practical stabilization and associated efficiency trade-offs may be 
drawn for more general target states. 

6 Outlook 

We have introduced a systematic method, following a system-theoretic approach, to determine whether 
a general mixed state of a finite-dimensional multi-partite quantum system may be the unique fixed 
point for a natural class of QL FF Markovian dynamics - for given locality constraints. If appropriate 
conditions are obeyed, a constructive procedure to synthesize stabilizing dynamics is provided. For 
states that are not full-rank, the same necessary conditions hold, however a method for constructing 
stabilizing generators is provided only if an additional “support” condition is satisfied. In both cases, 
stabilization may be achieved without requiring Hamiltonian control resources in principle. We have 
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presented a number of QIP and physically motivated examples demonstrating how our tools can nat¬ 
urally complement and genuinely extend available techniques for fixed-point convergence and stability 
analysis - including QL stabilization of Gibbs states of non-commuting Hamiltonians (albeit for finite 
system size). Altogether, beside filling a major gap in the existing pure-state QL stabilizability analy¬ 
sis, we believe that our results will have direct relevance to dissipative QIP and quantum engineering, 
notably, open-system quantum simulators. 

A number of research questions are prompted by the present analysis and call for further investi¬ 
gation. Determining whether our necessary condition for FFQLS is, as we conjecture, always sufficient 
on its own even for non-full-rank states is a first obvious issue to address. From a physical standpoint, 
in order to both understand the role of Hamiltonan control and to make contact with naturally oc¬ 
curring dissipative dynamics, it is important to scrutinize the extent to which the stabilizing dynamics 
obtained in our framework may be compatible with rigorous derivations of the QDS - in particular in 
the weak-coupling-limit, where the interplay between Hamiltonian and dissipative components is crucial 
and demands to be carefully accounted for mm Still within the present QL stabilization setting, an 
interesting mathematical question is to obtain a global characterization of the geometrical and topolog¬ 
ical properties of the FFQLS set, beginning from pure states. As alluded to in the last remark above, 
answering these questions may also have practical implications, in terms of approximate QL stabiliza¬ 
tion. Related to that, while our main focus has indeed been on exact asymptotic stabilization, it may 
also be beneficial (possibly necessary) to tailor approximate methods for analysis and/or synthesis to 
specific classes of states. A natural starting point could be provided here by graph states, which have 
recently been shown to arise as arbitrarily accurate approximations of ground states of two-body FF 
Hamiltonians |18) . 

From a more general perspective, as remarked in the Introduction, the present analysis fits within 
our broader program of understanding controlled open-quantum system dynamics subject to a resource 
constraint. In that respect, an important next step will be to tackle different kinds of constraints - 
including, for instance, less restrictive notions of quasi-locality, which allow for exponentially decaying 
interactions in space (as in [29)1: or possibly reformulating the QL constraint away from “real” space 
and the associated tensor-product-decomposition, but rather relative to a preferred operator subspace, 
in the spirit of “generalized entanglement” [4]. Lastly, and perhaps somewhat counter-intuitively, as 
our results on separable non-FFQLS evidence, it is interesting to acknowledge that the QL stabilization 
problem need not be trivial for classical probability distributions. For both the classical and the quantum 
setting, we find it tempting to speculate that interesting mathematical connections may exist with the 
quantum marginal problem and general quantum joinability EDllH]- 
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Appendix A: Equivalence between FFQLS and DQLS for pure 
fixed points 

Proposition A.l. (Pure-state FFQLS is DQLS) A pure state p = DQLS if and only if 

it is FFQLS. 


Proof: If |^/’)('0| is DQLS, then by definition there exists QL dynamics with no Hamiltonian component 
and QL noise operators satisfying Dk\if) = 0 in standard form. Then, clearly, both £(|^)(i/)|) = 0 and 
= 0 for each neighborhood, making |^) FFQLS. 

To prove the converse implication, note that it follows from Theorem [S] (proved in [^) that a pure 
state |'0)('!/'| is DQLS if and only if span(|'!/))('0|) = H(SjVj (I?/’)) We show that this equality 

follows if |V^)('0l is FFQLS. If so, then Theorem [7] implies that 

span(|V')(V'|) = 

3 

The Schmidt span of a pure density operator |■^/')('0I is simply Etv'j(I'*/’)('*/'I) = 'S[supp(pjVj)]■ Then, 
H[supp(pjv'j)] is already a modular-invariant -distorted algebra, thus we do not need to enlarge it 
to obtain (^At^ (|■*/')(l/’l)) With this, and the fact that B{V) ®B{W) = B{V ® W), we 

may write 

span(|V’)(V'|) = f]B{swpi>{p^r^ ® %,))• 

3 

Lastly, we use the properties that Cl^BfVi) = B{f].^Vi) and B{V) = B{W) V = W to reduce the 
above equality from operator space to vector space, 

supp(i'0)('i/'i)= 

3 

implying that the target state Itf) is DQLS, as desired. □ 


This result yields a characterization of FF dynamics which uniquely stabilize a pure state: 

Corollary A. 2. If a pure state p = \ip){'ip\ is FFQLS, then there exists a generator that makes |'0)('!/'| 
FFQLS and is of purely-dissipative form: 

Ci-) = J2Dk-Di-^{DiDk,-}, 

k 

where each Lindblad operator Dk is a neighborhood operator and, in addition, obeys Dk\il3) = 0 for each 
k. 


Appendix B: A result on fixed states of convex combinations of 
CPTP maps 

Theorem B.l. (Common fixed points of sums of CPTP maps) Let T = 'YhkPklk be a sum of 

CPTP maps, with pk > 0 and ^j^Pk = 1; cind assume that both the following conditions hold: 

(i) there exists a common fixed point p £ fix(7fc) for all k; 
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(ii) supp(p) = supp(ker(T) = supp(fix(7fc)) for each k. 

Then p' is invariant under T only if it is invariant under all 7k, that is: 


p' e fix(T) p' G fix(7^). 


Proof: We show that fix(T) = P|j.fix(7fc). First, linearity of T implies that fix(T) > P|j,fix(7fe). The 
hypotheses are needed to show fix(£) < P|^fix(7fe). By (ii), p is a maximal-rank state in fix(T) and 
p G fix(7fe) for all k. Thus, its support 71 is invariant for all maps and we can construct reduced maps 
Tfe as we did in Theorem O for reduced generators C. The reduction of p to this support, call it p, is 
full-rank. In view of (ii). Theorem [T] implies that 

fix(T) = p5(fix(T^) 0 0)p5 and fix(7fc) = (fix(7)J) 0 0)p^, 

for all k. Then, by Theorem [I] we also have that 

p5fix(Tt)pl = p^alg{T}V^ and p^fix('7^^)p^ = p^alg{7fe}V^- 

Since alg{T} > alg{7fc} for all k, the relevant commutants satisfy alg{T}' < alg{7fc}', for all k. This 
inequality may be used to bridge the previous equalities, yielding 

fix(T) = p5 (fix(T'^) 0 0)p^ = p5alg{T}'p^ ® O 

lA 

fix(7fc) = p5(fix('^''')0O)p^ = p5alg{7fe}'p^ ® O, 
for all k. From this we obtain fix(T) < Hfc fix(7fe), which completes the proof. □ 

Appendix C: Results and proofs for pure Dicke states 

Proposition C.l. (Pseudo-pure Dicke states) Pseudo-pure versions of Dicke states are not FFQLS 
in general. 

Proof: Consider, in particular, the following state of the general form given in Eq. (1551) : 

Ppp = (l-e)|(OOll))((OOll)|0eI/2^ 

To calculate the operator Schmidt span of this state, we first calculate that of the pure state. The 
l|234-Schmidt decomposition of |(0011)) is 

1(0011)) = ^(|O)|(O11))0|1)|(OO1))), 

with the analogous 123|4-Schmidt decomposition being obtained by symmetry. With the above pure- 
state Schmidt decomposition, we easily calculate that of the pseudo-pure state: 

Ppp = I/V2 0[e|(OOl))((OOl)|0e|(Oll))((Oll)|0(l-e)I/8] 

0eZ/^0 (|(011))((011)| - |(001))((001)|)/V2 

+eY/y/2 0 (i|(001))((011)| - t|(011))((001)|)/y2 
+eX/V2 0 (|(001))((011)| 0 |(011))((001)|)/V2 
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Our prescription is to distort the 234-Schmidt span by /9234 calculate the generated *-algebra. The 
distortion takes the four Schmidt basis operators above to (after rescaling) 

1234 , |( 011 ))(( 011 )|-|( 001 ))(( 001 )|, 

*|( 001 ))(( 011 )| - *|( 011 ))(( 001 )|, |( 001 ))(( 011 )| + |( 011 ))(( 001 )|. 

The *-algebra generated by the span of these operators is , where 'Hw = span{|(001)), |(011))} 

and is its orthogonal complement. Re-distorting this *-algebra with P2341 that it is un¬ 
changed by such an action. Further, it is invariant under Af 1. Appending the minimal fixed- 

point set corresponding to the 1|234 bipartition is 

® .7^p234 (^234(Ppp)) = Bi^i ® 'Hw) ® {H{Hi) ® ^H^)- 

By symmetry, the minimal fixed-point set corresponding to the 123|4 bipartition is 

■^Pi23(^123(Ppp)) ® B{H 4 ) = B{Hw ® H 4 ) ® ® B{H 4 ))- 

While the two minimal fixed-point sets both contain the target state, they also both contain the identity 
operator. Their intersection is thus the two-dimensional space 

B[Hi ) ® ■^P234 (^234(Ppp)) ■Api23 (^123(Ppp)) ® B[H 4 ) = span{I, Ppp}, 
violating the necessary condition for FFQLS. □ 

We next present the proof of Proposition [T] The following two lemmas are useful to characterize 
and constrain the intersection of Schmidt spans for generalized Dicke states. 

Lemma C.2. Consider a pure state \tjj), neighborhood structure {Mj}, with j = 1,... ,L, and corre¬ 
sponding Schmidt spans STv^d^/)). For each neighborhood, let Gj be a representation of a group which 
acts trivially on EjVj (|4')). For each such group, we trivially extend its action from H^f^^ to Hjg'^ ® 
by tensor-factoring to its elements. Then, the group generated by all the neighborhood groups, 

G= (Gi,G2 ,...,Gl), 

must act trivially on the intersection 

Proof: Let v £ T,jy. {\tj})) 0 Hjj'.. Then v £ Ejy. {\tfj)) 0 Hj^. for all j = 1,..., L. It then follows, by 
assumption, that the elements of each Gj act trivially on v. □ 

Lemma C.3. For any generalized Dicke state and a connected neighborhood structure, the intersection 
of all neighborhood Schmidt spans is contained in the symmetric subspace. 

Proof: By Eq. (|48l) . it follows that any Schmidt span of a generalized Dicke state is spanned by 
Dicke states, and hence is contained in the symmetric subspace. By definition, the latter is the trivial 
representation of the symmetric group. Thus, for any Dicke state, the Schmidt span with respect to 
some neighborhood is acted on trivially by permutations of that neighborhood’s subsystems. For a 
connected neighborhood structure, the set of neighborhood-wise permutations generates the full 
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By Lemma C.2, the intersection of Schmidt spans for that neighborhood structure must be acted on 
trivially by whereby the claim follows. □ 

Proposition V.8. (FFQLS Dicke states) Given n qudits and a connected neighborhood structure 
M, there exists a (non-factorized) FFQLS generalized Dicke state relative to M if d{m — 1) > n, where 
m is the size of the largest neighborhood in Af. 

Proof: We show that the entangled n-qudit Dicke state Dn,m defined in Eq. (|4^ in the main text is 
FFQLS relative to any connected neighborhood structure containing a neighborhood of size (at least) 
m. As established in Lemma C.3, the assumed connectivity of Af implies that the intersection of the 
extended Schmidt spans 'Hjf lies in the symmetric subspace. Thus, it remains to show 

that any symmetric state not in span(D„_m) will fail to be included in some extended Schmidt span. 
Specihcally, we show that any such state fails to lie in the Schmidt span corresponding to a size-m 
neighborhood, which we label A/q. 

The candidate FFQLS Dicke state Dn^m has the special property that it contains no more than 
m — 1 of symbols 0,... ,d — 1, and no more than r of the symbol d — 1. Of all states in the symmetric 
subspace of n qudits, Dn^m is the only one satisfying this property; viewing each label 0 ,..., d — 1 as 
a bosonic mode, Dn,m constitutes the only possible n-particle filling of these modes without exceeding 
the imposed occupancy limits of m — 1 (and r for mode d — 1, respectively). It remains to show that 
any symmetric state vector \4i) not satisfying this mode-occupancy-limit property, and hence not in 
span(D„_m), fails to be in at least one extended Schmidt span of the state Dn^m- A symmetric state 
which fails to satisfy the above mode occupancy limits will admit an orthogonal decomposition of the 
form 

\(f) = a| ® +/3k), 

m 

if a mode j = 0 ,..., d — 1 is over-occupied, or of the form 

|r) = 7| d - 1, d - 1 ,... d - 1 .. .)7yo 0 +<5|i), 

r+1 

if mode d — 1 is over-occupied. The extended Schmidt span 'Fj^g{Dn,m) 0 'fLjJg i® spanned by vectors 
of the form |A)7yg 0 \i)jj., where |A)7y(, are Dicke states with m symbols drawn (without replacement) 
from the symbols of Dn^rm and the |z) form a basis for Td-jj.. Crucially, Dn^m is such that no |A) will 
have more than m — 1 occurrences of j, implying that 

{Mjj =0: (A| d- l,d - 1,.. .d- 1.. .)aAo = 0. 

r+1 

From this fact, the length of both \(j>) and |r) will be decreased when orthogonally projected onto the 
extended Schmidt span (more precisely, they cannot exceed the lengths 11/31^)11 < 1 and ||d|t)|| < 1, 
respectively). Hence, the states |(/)) and |r) are not in the extended Schmidt span A/q. This shows that 
F>n,m is the only state vector in the intersection of its extended Schmidt spans, and hence it is FFQLS, 
as claimed. □ 
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